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We investigate theoretically the effects of intralayer and 
interlayer exchange in biased double-layer electron and hole 
systems, in the absence of a magnetic field. We use a varia- 
tional Hartree-Fock-like approximation to analyze the effects 
of layer separation, layer density, tunneling, and applied gate 
voltages on the layer densities and on interlayer phase co- 
herence. In agreement with earlier work, we find that for 
very low layer separations and layer densities, an interlayer- 
correlated ground state possessing spontaneous interlayer co- 
herence (SILC) is obtained, even in the absence of interlayer 
tunneling. In contrast to earlier work, we find that as a 
function of total density, there exist four, rather than three, 
distinct noncrystalline phases for balanced double-layer sys- 
tems without interlayer tunneling. The newly identified phase 
exists for a narrow range of densities and has three compo- 
nents and slightly unequal layer densities, with one layer being 
spin-polarized and the other unpolarized. An additional two- 
component phase is also possible in the presence of sufficiently 
strong bias or tunneling. The lowest-density SILC phase is the 
fully spin- and pseudospin-polarized "one-component" phase 
discussed by Zheng and co-workers [Phys. Rev. B 55, 4506 
(1997)]. We argue that this phase will produce a finite in- 
terlayer Coulomb drag at zero temperature due to the SILC. 
We calculate the particle densities in each layer as a function 
of the gate voltage and total particle density, and find that 
interlayer exchange can reduce or prevent abrupt transfers of 
charge between the two layers. We also calculate the effect of 
interlayer exchange on the interlayer capacitance. 

PACS numbers: 71.45. Gm, 73.20.Dx, 73.40. Gm, 73.40.Kp 



I. INTRODUCTION 

In the last several years, double-layer electron and hole 
systems have provided an exceptionally useful tool for in- 
vestigating the effects of interparticle Coulomb interac- 
tions in two dimensions, particularly at low particle den- 
sities where exchange and correlation effects are signifi- 
cant. Tjius has been especially true in the quantum Hall 
regime,Era where the combination of a strong perpendic- 
ular magnetic field (which quenches the kinetic energy) 
and very small layer separation (which enhances inter- 
layer exchange) stabilises remarkable interlayer-coherent 
quantum Hall states.uEi Even in the absence of any mag- 
netic field, interlayer capacitance and drag measurements 
in double-layer electron systems (2LES's) have provided 
quantitative measures of the effects of— electronic inter- 
actions on both, the thermodynamics^ and transport 
(Coulomb dragB) of two-dimensional electron and hole 



systems. Unless otherwise specified, we shall take the 
notation 2LES to also include double-layer hole systems. 

Our work is motivated by the rapid pace of advance- 
ments in the engineering of double-layer semiconductor 
devices. We expect that high-mobility double-layer de- 
vices will eventually be built with both (1) separately 
contactable layers, and (2) layer separations and carrier 
densities so small that the interlayer correlations between 
the carriers are substantial, perhaps even without the aid 
of a strong quantizing magnetic field. Such devices will 
allow direct measurements of the effects interlayer many- 
body effects in double-layer systems. As a starting point 
to analyze the zero magnetic field situation, we have de- 
veloped a simple mean-field model that incorporates both 
intralayer and interlayer exchange in biased double-layer 
electron and hole systems in the absence of a magnetic 
field. We use the model to calculate theoretically the 
effect of interlayer exchange on the layer densities and 
interlayer capacitance as a function of layer spacing, par- 
ticle density, and applied gate voltage. 

Gated double-layer systems have the great advantage 
of allowing the layer densities of electrons or holes to 
be varied by the application of a bias (gate voltage). At 
high densities, the kinetic energy per unit area dominates 
the exchange and correlation energies, and in a transla- 
tionally invariant system, the symmetry and properties 
of the ground state of a many-body system are in one- 
to-one correspondence with those of a free-electron gas. 
At lower electron densities, the Coulombic exchange and 
correlation energies can produce qualitative changes in 
the nature of the many-particle ground state: numerical 
work on two- and three-dimensional electron gases show 
that a spin ferromagnetic state is obtained at low densi- 
ties, which is eventually supplanted by a Wigner crystal 
state at the lowest densities. The low-density ferromag- 
netic state of the interacting electron gas was anticipated 
some 70 years ago by Blochjlj ferromagnetism can be 
found even within the Hartree-Fock approximation when 
the exchange interaction energy, which favors occupation 
of single-particle states of the same spin, dominates the 
kinetic energy, which favors reducing the Fermi energy 
by equal occupation of both spin states. 

Multilayer semiconductor devices enhance the effects 
of interparticle interactions through the combination of 
reduced dimensionality and low particle density, and by 
the presence of an additional electronic degree of free- 
dom, the layer index.Ej In double-layer systems, layer oc- 
cupancy can be specified by a introducing a pseudospin 
variable thai points up for one layer and down for the 
other layer S3 Extending the notion of an exchange-driven 
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ferromagnetic transition to double-layer systems suggests 
that at low enough densities, the electronic ground state 
should be both spin and pseudospin polarized. Such 
reasoning, supported by Hartree-Fock calculations, led 
Ruden and Wu to propose that, at low enough electron 
densities and for small enough layer separations, the elec- 
trons of a 2LES would minimize their ground-state en- 
ergy by having all electrons occupy a single layer Jl3 The 
Ruden- Wu scenario implies that as the total density of a 
balanced 2LES is lowered, the electrons should eventually 
experience an interlayer charge transfer instability that 
spontaneously empties out one of the two layers. Ruden 
and Wu also suggested that the inherent bistability of the 
resulting low-density 2LES (either layer could be the one 
to lose or gain particles) would constitute an exchange- 
driven logic gate. The possibility of an exchange-driven 
interlayer charge- tpinsfer instability has fc>3§fti a subject 
of both theoreticalEJta and experimentalll3~E3 interest. 

Although pseudospin polarization at sufficiently low 
densities is likely, it does not require unequal layer den- 
sities. This has been demonstrated theoretically in great 
detail for closely spaced double-layer systems in a strong 
magnetic field at unit filling factor and appears to give 
a good explanation of experimental results-EHa The key 
point is that electrons in double-layer structures are not 
restricted to occupying only one of two layer eigenstates: 
quantum mechanics allows states that are superpositions 
of the two layers. The layer pseudospin must there- 
fore be treated as a Heisenberg variable, as was done 
in Refs. || and in the quantum Hall regimer-iand by 
Zheng and co-workers in zero magnetic fieldJiH rather 
than as an Ising variable (where only "up" and "down" 
allowed) as was done by Ruden and Wu. For example, 
when interlayer tunneling is present, the single-particle 
eigenstates are symmetric and antisymmetric combina- 
tions of layer states. A major insight of Refs. |^-^ 
was the concept of "spontaneous interlayer coherence" 
(SILC): at sufficiently small layer separations, electrons 
can spontaneously create and occupy linear combinations 
of layer states in which each layer has the same average 
number of electrons, even without any interlayer tun- 
neling. The spontaneous formation of superposed layer 
states in double-layer quantum Hall (2LQH) sytems can 
be accomplished by the interlayer exchange interaction 
alone. SILC in balanced 2LQH systems corresponds to 
XY pseudospin ferromagnetism in which the pseudospins 
spontaneously magnetize, but do not point either up or 
down, since neither layer has (on average) more particles 
than the other.H 

The application of SILC to the zero magnetic field 
case was first made by Zheng and co-workers, til who pQn- 
sidered the same model system as Ruden and WuE3 - 
electrostatically balanced zero-thickness layers of inter- 
acting electrons without interlayer tunneling - but came 
to a very different conclusion. They proved within the 
Hartree-Fock approximation (HFA) that (1) the 2LES 
becomes spin ferromagnetic before it becomes pseudospin 
ferromagnetic for any finite layer separation (this possi- 



bility was not considered by Ruden and Wu), and (2) at 
low enough densities and small enough layer separations, 
the pseudospin ferromagnetic state possesses SILC, with 
all electrons occupying one subband composed of a su- 
perposition of layer states with equal average density in 
each layer. Conti and Senatore have also argued that for 
electrostatically balanced layers, the 2LES ground state 
cannot be one in which all electrons are eigenstates of 
the same single layenl3. SILC at sufficiently small layer 
separations should also follow from earlier calculations 
presented in Ref. [2^, although no explicit inference of 
SILC was made in that work. Recent work that goes be- 
yond the HFA and includes correlation effects within the 
STLS approximationcj also finds that SILC is favored 
over single-layer occupancy for balanced layers that are 
sufficiently close together. c3 

When studying a 2LES using a density-functional ap- 
proach, it is important to note that SILC is a nonlo- 
cal effect. Calculations based on local-density approxi- 
mations will not find SILC if they treat the pseudospin 
as an Ising-like variable O The same caveat applies to 
the work of Ruden and Wu, who used a restricted HFA 
that excluded the possibility of SILC. Once the intralayer 
separation between electrons becomes comparable to the 
interlayer separation between layers, the possibility of 
interlayer correlations (such as SILC) must be consid- 
ered. The spontaneous charge-transfer state predicted 
by Ruden and Wu follows quite generally (even beyond 
the HFA) from the fact that when interlayer correlations 
are ignored, the negative compressibility of the electron 
gas guarantees an interlayer charge-transfer instability 
at sufficiently small layer spacing. But it is precisely at 
small layer spacings that interlayer correlations become 
important, and so their effects must be included to obtain 
physically meaningful results. 

We have extended the study the effects of Coulom- 
bic exchange in double-layer electron and hole systems 
to include an applied bias due to front and back gate 
voltages, while allowing for the possibility interlayer ex- 
change. In the balanced case, we have found that the 
four- to two-component transition is always interrupted 
by the presence of a three-component phase with slightly 
unequal layer densities. There are therefore four possi- 
ble noncrystalline phases for a 2LES with balanced gates. 
Under bias or tunneling a second (pseudospin-polarized) 
two-component state is also possible. We have enumer- 
ated the transitions between the five allowed states in the 
presence of bias, and explored the effects of bias on the 
one-component state. 

The rest of this paper is organized as follows: In Sec. ||, 
we introduce a simplified model for double-layer systems, 
review the concept of interlayer capacitance, and give a 
general criterion for stability against spontaneous inter- 
layer charge transfer. In Sec. Ill, we develop a mean-field 
approximation for biased double-layer systems that al- 
lows for the possibility of interlayer coherence. In Sec. IV, 
we examine the special case of electrostatically balanced 
layers, enumerate the resulting four possible noncrys- 
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talline phases, and explore the onset of interlayer coher- 
ence and its effect on the size of the subband splitting. 
We also obtain a phase diagram for the balanced case, 
and perform an alternate calculation for the onset of the 
one-component phase. In Sec. [v|, we explore the effect of 
bias on the layer occupancies and interlayer capacitance 
for large, intermediate, and small layer separations. We 
develop simple models capable of closely fitting experi- 
mental layer-occupancy data, and expore the transitions 
between different phases induced by layer imbalance. In 
Sec. VI, we analyze the onset and properties of the one- 
component state under bias. We summarize our findings 
and speculate on the possible relevance of the se results 
to the strong magnetic-field regime in Sec. VII. 



II. 



DOUBLE-LAYER MODEL AND INTERLAYER 
CAPACITANCE 



In this section, we introduce an idealized model of 
double-layer systems. We review the condition for ther- 
modynamic equilibrium between the inner layers, ob- 
tain a necessary condition for stability against inter- 
layer charge transfer, and review an experimentally use- 
ful measure of the interlayer capacitance, the Eisenstein 
ratio£3 

Figure [j] illustrates schematically the geometry of the 
2LES device. We treat the quantum wells as zero- 
thickness layers sandwiched between two plates of neu- 
tralizing charge, which represent the effects of the front 
and back gates. The distance between the front gate (at 
far left) and the first layer (layer 1) is Dp', that between 
the back gate (at far right) and the second layer (layer 
2) is Db] the interlayer separation is denoted by d. Typ- 
ically, d ~ 10 nm, Dp ~ 1/im, and Db ~ 1 mm, so that 
d <C Dp <C Db\ thus, Fig. [l] is not at all to scale. 
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FIG. 1. Schematic figure of a double-layer device with 
front-gate surface charge density epF at left, areal charge den- 
sities —eni and — em in the quantum wells, and back-gate 
surface charge density eps at right. 

The two inner layers are assumed to be in thermody- 
namic equilibrium with each other, and the voltage of the 
front (back) gate relative to the common chemical poten- 
tial of the inner layers is denoted by Vf (Vb). It is also 
assumed that small changes SV a in the gate voltages V a 



(a = F, B) produce small changes in surface charge den- 
sities only at the gates (eSp a ) and in the layers (—eSrii, 
i = 1,2). Overall charge neutrality requires that the total 
charge density vanish: 



epp + eps — eni — eri2 = 0. 



(1) 



(Strictly speaking, we only require that the change in the 
total charge density vanish: SepF + SepB — 5en\ — 5eri2 = 
0.) In writing Eq. ([!]), we have assumed that any stray 
charge not included in Eq. ([!]) is unchanged when the 
gate voltage are varied. This implies that the only sig- 
nificant effect of the stray charges is to shift the gate 
voltages by constant (empirically determined) amounts. 
Sheet charge densities on the gates and inner layers pro- 
duce electric fields between the double-layer and the gates 
(E a ) and between the two layers (-E12) according to 
Gauss's law, 



E a = (e/e)p a , 
E12 = (e/e)(pF 



(2) 



7ii ) = (e/e)(n 2 ~Pb)- 



We now obtain the conditions for thermodynamic equi- 
librium and stability between the layers. Regarding the 
gate charge densities p a as fixed quantities, we seek the 
values of the particle density in the inner layers that min- 
imize the total energy per unit area £o/L x L y . Following 
Ref. ^6|, we separate the electrostatic part of the total 
energy per unit area from the rest: 



£0 



-E 12 d + e(n 1 ,n 2 ) 



(3) 



up to an irrelevant constant, where the first term is the 
interlayer electrostatic energy density, with E12 given by 
Eq. (|2|). The quantity e (^1,^.2) represents the total en- 
ergy per unit area for a fully interacting double-layer sys- 
tem in which each layer contains a uniform neutralizing 
charge density, for particle densities n\ and n 2 in layers 1 
and 2, respectively. The Fermi, exchange, and correlation 
energies, both intralaycr and interlayer, are contained in 
e(n 1 ,7i 2 ), and it is thiS|-quantity that is calculated using 
many-body techniques .c3 In the next section, we make an 
approximate calculation of £(77,1,712) which includes the 
effects of interlayer and intralayer exchange. 

To obtain the condition for thermodynamic equilib- 
rium between the inner layers, we note that for fixed ex- 
ternal charge densities p a , the constraint of overall charge 
neutrality implies that the particle density in one layer 
(e.g, Tii) is determined by that in the other layer (e.g., 
7i 2 ): 77.1 = Pf + Pb — 7i2- We may thus regard the total 
energy per unit area £q/L x L v as a function of the layer 
density 77-2, and extremize j^/j^Ly with respect to 7i 2 at 



fixed pf and ps to obtainLil'E 



where 



Mi - M2 = eE 12 d, 



fii = ds(ni,n 2 )/dn i 



(4) 



(5) 
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is the chemical potential measured relative to the energy 
minimum of layer i. Equation (ji|) states that the dif- 
ference in the layer values of the chemical potential is 
equal to the electrostatic potential energy difference be- 
tween the layers. If the equation of state determining 
fii(rii, ria) were known, then Eqs. @) and (Q) would to- 
gether determine the values of layer densities n\ and n 2 
for which the total energy £q/L x L v is an extremum. 

We now examine a necessary condition for interlayer 
thermodynamic stability (i.e., for the local extremum 
to be a local minimum of the energy). First we follow 
Rcf. [h] and introduce a set of lengths that describe the 
dependence of the layer chemical potentials fii on the 
layer densities nj, 



d 2 e 



e 2 drij e 2 dn^drii 



(6) 



For the extremum condition in Eq. (Q) to represent a local 
minimum of the total energy per unit area, we require 
that the second derivative of £q/L x L v with respect to n 2 
be positive. (We again regard n\ as being determined 
by n 2 for fixed p a by the requirement for overall charge 
neutrality.) This gives a necessary condition for stability: 



d + si + s 2 > 0, 



where 



s 2 = s 22 - s 2 i. 



(7) 



(8) 



The above inequality guarantees that the 2LES is at least 
metastable. If Eq. (^) is violated, then the layer densities 
constitute a local energetic maximum rather than a local 
minimum, and there will be an interlayer charge insta- 
bility that causes charge to flow between the layers until 
a new energetic minimum is reached. 

In the absence of interlayer correlations (the case con- 
sidered in Refs. [D] and |l^), s 12 = s 2 \ = 0, and the length 
Su is directly related to, the electronic compressibility Ki 
in layer i according too 



(9) 



Equations (]?]) and (^|) together with the experimentally 
measured negative compressibility (k < 0) of the elec- 
trons imply that when interlayer correlations are ignored 
(sx2 = 0), sufficiently small interlayer separation d will 
always lead to a charge-transfer instability when the den- 
sity of either (or both) layer is sufficiently small. This 
result is true even when intralayer correlations are in- 
cluded beyond the HFA calculation of Ruden and Wu; 
only interlayer correlations that produce sufficiently neg- 
ative values of s\ 2 to satisfy Eq. ([?]) can suppress inter- 
layer charge-transfer instabilities at very small interlayer 
separations. 

If the two layers of the 2LES sample can be contacted 
separately, then the Eisenstein ratio Re provides a sensi- 
tive measure of the interlayer capacitance that avoids the 



large gate-distance faptors that dominate the gate capac- 
itances per unit area.E3c3 The Eisenstein ratio is defined 
as the ratio of the differential change in the electric field 
E\ 2 between the inner layers to that of the electric field 
Ep between the front gate and the inner layers: 



Re = 



SE 



12 



5E f 



1 — Sni/dp 



F, 



(10) 



where we have made use of Gauss's law, Eq. (g). In the 
following sections, we calculate the layer occupancies rii 
as a function of the gate charges p a ; we then use Eq. jlfj| ) 
to obtain Re by computing the derivative of n\ with 
respect to pp. In the classical limit (corresponding to 
large enough particle densities and layer separations so 
that only the electrostatic energies are relevant), n\ = pf 
(for pf > 0), so that by Eq. @, R E = 0. Note also that 
if m = (e.g., due to pp < 0), then Re — 1- By using 
Eqs. (^) and (|J) to express differential changes in Ei 2 in 
terms of the electronic lengths sy and using Gauss's law, 
the Eisenstein ratio may be expressed as 



Re = 



si — s 2 5Eb /SEp 



d + si 



•S-2 



(11) 



The Eisenstein ratio has an especially simple form 
for fixed total density since then 6Eb — —5Ep; from 
Eq. (|l|), Re = s/(d+ s), where s = si + s 2 . However, 
most experiments fix the back-gate voltage Vb rather 
than the total density, and sweep the front-gate voltage 
Vf- Because of the relatively large size of the back-gate 
distance Db, the constraint of fixed back-gate voltage 
(SVb — 0) is very nearly equivalent to that of fixed back- 
gate sheet charge density {&pb = or, by Gauss's law, 
8E B = 0). For SE B = 0, Eq. @ shows that the Eisen- 
stein ratio is very nearly 



R E w (6E 12 /SE f ) % 



si 



si + s 2 



(12) 



The advantage of measuring Re rather than the usual 
gate capacitances per unit area (such as eSpF / SVf) is 
that Re depends only on the electronic lengths Si and 
the interlayer distance d, not on the much larger gate 
distances D a . The difficulty in measuring Re is that 
(at least one of) layers 1 and 2 must be separately con- 
tactable, which becomes increasingly difficult as the layer 
separation d becomes very small. Nonetheless, measure- 
ments of the Eisenstein ratio have been used to demop. 
strate the negative compressibility of the electron gas,E2l 
and it is to be expected that devices with separately 
contactable layers will be built with increasingly narrow 
layer separations. Note also that Re is very sensitive 
to charge-transfer instabilities. In fact, from Eq. (0), 
the condition for the onset of an instability to interlayer 
charge transfer, d + s± + s 2 = 0, shows that Re formally 
diverges at the instability. This can also be seen from 
the relation Re = 1 — Sni/Spp, since a charge transfer 
instability would produce an abrupt change in the layer 
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density n\ in response to a small change in the gate den- 
sity pf- Although evidence for abrupt interlayer charge 
transfers has been renorted based on Shubnikov-de Haas 
(SdH) measurements JL3 the Eisenstein ratio would be a 
far more sensitive measure of abrupt interlayer charge 
transfers. 

In subsequent sections, we calculate the layer densities 
Tii as a function of the front-gate particle density pf, for 
fixed back-gate particle density pb- Fortunately, pb may 
be found experimentally from SdH measurements as the 
value of the layer densities when the system is balanced: 
i.e., for equal layer densities (n\ = n%) and minimum sub- 
band separation (n a — rib). Once ps (which we assume 
is very nearly constant) is known, pp may be determined 
from charge conservation by measuring the total layer 
density, either by SdH (as the sum of the subband densi- 
ties) or by Hall effect measurements. It is therefore possi- 
ble to determine pf experimentally, without recourse to 
the gate voltages. Of course, experimentally, it is Vp that 
is varied directly while Vb is kept fixed, and pp changes 
in response to Vp (while ps changes very little for large 
Db)- We describe how the gate voltages may be deter- 
mined from a knowledge of pf,Pb, and fj,i in Sec. A 1 of 
the Appendix. 

III. MEAN-FIELD APPROXIMATION 

In this section, we introduce the microscopic Hamilto- 
nian for the gated double-layer system and make a varia- 
tional approximation for the ground-state wave function 
that allows for interlayer coherence. The resulting ap- 
proximate ground-state energy per unit area depends on 
both intralayer and interlayer exchange. 

We now consider the microscopic Hamiltonian for the 
double-layer system illustrated in Fig. |l|. We idealize the 
inner layers as being two-dimensional, and treat inter- 
layer tunneling in the tight-binding approximation. The 
Hamiltonian for the interacting system is then 

H = E £kC ]*s C 3^ - i E^lks^ks + 4ks C lks) (13) 



j'ks 



ks 



2rjr y Y. E E Yhhb) 

q jlkisi J 2 k2S 2 



XC. 



iki+qsi C j2k 2 -qs2 C J 2k2S2C J lklSl 



~ E E V 3°(9 = °K C lk S C ik S 
jks a. 

+ E V <*M = °)P°>Ppi 

where Cyks ( c jks) denotes the second-quantized destruc- 
tion (creation) operator for an electron or hole in layer 
j with momentum Kk and spin s. Here = Ti 2 k 2 /2m* 
is the kinetic energy in the effective mass (to*) approxi- 
mation, t is the interlayer tunneling amplitude, and the 
Fourier-transformed Coulomb potential is given by 



2eq 



-qdij 



(14) 



where dij is the distance between layer i and layer j. The 
indices a, (3 in Eq. ( |l4| ) take the values F (front gate) 
or B (back gate), V a p is the direct Coulomb interaction 
between gates a and /?, and V a j is the direct Coulomb 
interaction between gate a and layer j. The last two 
terms in Eq. (Q) represent the direct interaction between 
the layers and gates, and between front and back gates, 
respectively. 

We use a mean-field approximation (MFA) that is vari- 
ationally based and that reduces to the Hartree-Fock ap- 
proximation for the balanced case of equal layer densities. 
This approximation includes the effects of interlayer cor- 
relations in the simplest possible way. The variational 
ground-state wave function is composed of two subbands, 
a and b, containing spin up (f) and spin down Q) elec- 
trons: 

/c4<fcfa| k3<kbi ^2^k a i fci</e a ^ 

i*o>= n' & L n b L n «l n <ti°>. 

k 4 k 3 k 2 ki 

(15) 

where k as and kb s denote the Fermi wave vectors for elec- 
trons or holes of spin s in subbands a and b. The creation 
operators for the subbands are related to the layer cre- 
ation operators by a canonical transformation that we 
take to be of the form 



"ks 



= cos(6/2)c\ ks + sin(0/2)e^ ( 



.t 

2ks 



(16) 



- sin(6»/2)e-^c t lk5 . + cos(0/2)c| 



:2 k. - 



When 9 — -k/2 and cf> = 0, subband a is the symmet- 
ric subband and subband b is the antisymmetric sub- 
band. In the language of pseudospin, the superposition 
of layer states in Eq. ( |l7| ) corresponds to treating the 
layer pseudospin as a Heisenberg, rather than an Ising, 
spin variableO The form of the canonical transforma- 
tion in Eq. (|l7|) is not completely equivalent to a fully 
self-consistent Hartree-Fock calculation because we have 
taken 9 and 4> to be independent of the wave vector k 
and spin s. It would be interesting to explore the ef- 
fects of including the k and s dependence of 9 and (f> in 
a future calculation. Our simpler variational calculation, 
which is equivalent to Ref. 17 for the special case of bal- 



anced layers, offers a reasonable starting point, which is 
probably qualitatively correct over a large range of layer 
densities. It certainly gives layer densities that are in 
close agreement with experimental values obtained from 
SdH measurements, as we shall show. 

The layer occupation numbers may be expressed in 
terms of the subband occupation numbers by using 
Eq. @: 

cos 2 (0/2)(4 s a ks ) +sin 2 (0/2)(&j cs 6 ks ) (17) 



( c lks c lks/ ^/-"A^ka"**/ T^iu V/^\"ks 

(4sC2ks) = sin 2 (0/2)(4 s a ks ) + cos 2 (0/2) (^^k,) 



5 



C2ks 



(4 ks ci ks ) = (4 ks c 2ks ) 



sin(0/2)cos(0/2)e i *((4 s a ks ) - (bt s b ks )) 



where the asterisk denotes complex conjugation, and we 
have used Eq. ( fL5| ) to eliminate cross terms such as 
(a ks &ks)- SILC occurs when 



q Us c 2 k S ) ^ o 



(18) 



in the absence of interlayer tunneling. According to Eqs. 
( |l8| ) and ([jJI), SILC requires that the following occurs: 
(1) 9 7^ 0, 7r so that the subband densities are different 
from the layer densities. SILC is therefore excluded when 
the pseudospin is treated as an Ising variable. (2) n a ^ n b 
so that the subband densities are not equal (nonzero 
pseudospin polarization). When n as = nj, s (completely 
unpolarized pseudospin), the MFA ground state can be 
expressed as the product of two uncorrelated single-layer 
wavefunctions by performing a global pseudospin rota- 
tion. 

Although it is straightforward to generalize our ap- 
proach to finite temperature, we calculate numerical re- 
sults in the limit of zero temperature (T = 0), both 
for the sake of simplicity and because measurements can 
(and have) been made on double-layer systems at low 
temperatures, even down to millikelvin temperatures in 
the quantum Hall regimefl For the zero-magnetic-field 
case treated here, we expect that finite temperature will 
not produce signifcant qualitative changes in the layer 
densities if the temperature T is below a fraction of the 
Fermi temperature Tp = Ep / ks , where Ep is Fermi en- 
ergy and ks is Boltzmann's constant. For a layer density 
of n = 10 10 cm -2 , Tp is roughly 4 K for n-type GaAs, 
and IK for p-type GaAs. The scale of the Hartree charge 
transfer energy, e 2 dn/2e, is larger than the Fermi energy 
except for ultrasmall layer separations (d < 5 nm for 
n-type GaAs and d < 1 nm for p-type GaAs.) 

The other reason we work at zero temperature is to ad- 
dress matters of principle, such as whether an interlayer 
charge transfer instability can occur when the layers are 
very close together; finite temperatures would presum- 
ably smear out such a transfer, if it could occur. In the 
limit of zero temperature, Eq. (|l5|) implies that 



(4 s a ks ) = e(k as - k), (bt s b^ s ) = e{k bs - k), (19) 

where the subband Fermi wave vectors and number den- 
sities are related through 



k as = y/4Tm as , k bs = 



(20) 



Summing Equations (|18j) over wave vector k relates the 
number densities of the layers to those of the subbands: 



nis + n 2s = n as + n bs , 

nis - n 2s = {n as - n bs ) cost 



(21) 



We may use the preceding equations to express the 
ground-state energy per unit area in terms of the subband 
occupancies n as and the angle 9: 



L T L U 



t(n a — rib) sin 9 cos <fi 



(22) 



— [(n a - rib) cos 9 — (pf — Pb) 
oe 

x y 



qas 



sin 2 9 s—^ 
AL T L lt ^ 



[Vu(q)-Vu(q)] 

x [Iaas(q) + hbs{q) - 2I abs (q)} , 



where vq — m*/(irfi 2 ) is the density of states per unit 
area for noninteracting spin-1/2 particles in two dimen- 
sions, q = ki — k2, and 



Ia(3s{q) 



1 



K 



Q(k as - 
x8(fc (3s 



|K + q/2|) 
-|K-q/2|). 



(23) 



Here K = (ki + k2)/2, and the subband indices a and 
j3 can be either a or b. Equation ( pi| ) says that I a p s (q) 
is l/(27r) 2 times the common area of two circles of radii 
k as and kp s whose centers are separated by q. When 
(3 — a, then kp s — k as , and the first exchange integral in 
Eq. (p3|) may be carried out explicitly: 



1 



2L X Ly 



^2Vn(q)I aas (q) = -- 



Aire 



,3/2 



(24) 



Equation (J2J) is just the exchange energy per unit area 
for a uniform single-layer spin-polarized two-dimensional 
electron gas of areal density n as . 

The last term in Eq. (^3|) , which contains the interlayer 
exchange contribution, may be conveniently expressed as 



sin 



'-(n a - n b ) 2 T, 



(25) 



4 2e 

where T is the interlayer exchange parameter, given by 



1 



L X L, 



E 

qs 



2 d/2e 



Iaas(q) + hbs(q) - 2l a bs(q) 



(26) 



n b ) 2 



The properties of T are described in Sec. A 4 of the Ap- 
pendix. The last term in Eq. (23) must in general be 
evaluated numerically, although it vanishes at d — or 
when n as — n bs . It also vanishes if 9 is or 7T, in which 
case subband a is the layer (1 or 2) with most particles, 
while subband b is the layer with the fewest particles. 
Ruden and Wu implicitly treated the layer pseudospin 
as an Ising-like variable with and ir as the only al- 
lowed values for 9; in their approximation, the last term 
in Eq. (^) vanishes, and the interlayer effects we shall 
discuss here do not appear. 



6 



For definiteness, we take n a > rib, n a-\ > n a i, and 
< 9 < it. Our procedure consists of finding the values 
of n as and 9 that minimize the expected energy per unit 
area, Eq. U23 ) . Within our variational approximation, we 
find (as inEZTthe HFA) that the spins in a given subband 
are always either completely polarized (ferromagnetic at 
sufficiently low subband densities) or completely unpo- 
larized (paramagnetic at higher densities). Real systems 
probably possess intermediate polarization for a range of 
low densitiesOo For finite t > 0, the ground-state en- 
ergy per unit area, Eq. (|23|), is minimized for = 0. 
In the absence of interlayer tunneling, the ground-state 
energy per unit area is independent of </>, provided that 
4> is constant; for convenience we set = 0. The layer 
densities n\ s and n 2s may be obtained from n as and 9 
via Eq. (p2[). We begin our calculations in the next sec- 
tion, by considering the case of electrostatically balanced 
gates. 

IV. BALANCED GATES 

In this section, we consider the case of electrostatically 
balanced gates (j>f = Pb), beginning with zero interlayer 
tunneling. This was-.the situation originally considered 
by Ruden and Wu,£3 and more recently in Refs. p7|Jl8| 
and ^3]. For balanced gates, our approximation is equiva™ 
lent to the unrestricted HFA of Zheng and co-workers Jl3 
and except for our analysis of the three-component phase, 
most of our results agree with theirs. 

A. Zero tunneling 

The balanced case raises an important question of prin- 
ciple: can exchange and correlation effects alone, unaided 
by applied gate biases and unhindered by interlayer tun- 
neling, ever produce a ground state in which the densities 
of the inner layers are not equal? Based on a restricted 
HFA (which did not allow for SILC) Ruden and Wu pro- 
posed that for small enough layer densities and layer sep- 
arations, the answer is yes. Zheng and co-workers argued 
recently that an unrestricted HFA (which allows for, but 
does not mandate, SILC) gives the opposite answer.EZl 
Wc find that, except for a small region in density which 
supports a three-component phase that has a slight layer 
imbalance, the layer densities are equal when the gates 
are balanced. 



1. Zero layer separation 

In order to classify the four noncrystalline phases that 
we find for gate-balanced double-layer systems, it is use- 
ful to begin with the idealized case of zero layer separa- 
tion (d — 0). For d — 0, the Hamiltonian is invariant 
under spin rotation, pseudospin rotation, and the inter- 
change of spin and pseudospin. It is the same as the 
Hamiltonian for a four-layer system of spinless fermions 



with zero separation between all the layers: layers, sub- 
bands, and spins become interchangeable labels for the 
four components. The d — double-layer system is there- 
fore equivalent to a single-layer two-dimensional system 
of fermions with CP(3) symmetry. At d = 0, the in- 
terlayer Hartree energy is zero and the interlayer (V12) 
and intralayer (Vn) Coulomb interactions are equal. As 
a consequence, the variational energy in Eq. ( ^3| ) is inde- 
pendent of 9 and (for t = 0) of when d = 0, 
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and the MFA is equivalent to the HFA. 

For generality, we first consider an A-componcnt sys- 
tem, where N is twice the number of layers: N = 4 for 
double-layer spin-1/2 systems. At d — 0, the MFA lacks 
intercomponent correlations; thus, the total energy of the 
system is just the sum of the individual energies e(n as ) 
of each component. We can investigate the distribution 
of component densities in the MFA ground state in an 
A-component system by taking all but two of compo- 
nent densities to be fixed, and minimizing the energy of 
the remaining two-component system. If we label the 
two components we seek to minimize as 1 and 2, then 
according to Sec. [n], the condition for stable equilibrium 
(local minimum of the total energy) is 



fJ-(rii) = fi(n 2 ), s(m) + s(n 2 ) > 0, 



(28) 



where fJ-(rij) = de{nj)/drij is the chemical potential 
relative to the minimum energy of component j and 
s(rij) = dfi(rij)/drij is inversely proportional to the com- 
pressibility of component j. When e(rij) is the sum of the 
kinetic and exchange energies as in Eq. (|27|), then the 
MFA energy is minimized only if (1) both layer densities 
are equal (rii = 712), or (2) one or both layers are empty. 
There are no intermediate possibilities in the MFA. This 
d = result for N = 2 gives the results found by Ruden 
and Wu: when ri\ + n 2 is sufficiently large, the compo- 
nent densities are equal; when ni+ii2 is sufficiently small, 
there is an exchange-driven intercomponent charge insta- 
bility that empties out one of the components. In the 
absence of intercomponent correlations, we expect that 
at low enough densities the component compressibilities 
will be negative and that therefore one of the components 
will empty out, even if intracomponent correlations are 
included. In the A-component d = MFA ground state, 
any pair of layers either has equal density or has at least 
one of the layers empty. 

There are therefore N possible MFA ground states in 
an A-component system at d = 0, characterized by the 
number of components p that have nonzero and equal 
densities. The remaining N — p components have zero 
density. Defining the dimensionless average interparticle 
spacing per component by r s = 1/ \Jir{nT /A)a§ for an 
A-component system, where a — Air eh 2 /m*e 2 is the 
effective Bohr radius, we may write the energy per unit 
area for the "p-component" MFA ground state as 
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\nrjpj 



3/2' 



(29) 



In the limit N — > oo, p becomes a continuous variable 
that minimizes the energy for p = p^ , where 



Pe 



N 



3tt 
4r, 



7ra 2 n T , 



(30) 



where we have assumed that r s > (37r/4); otherwise, 
Poo = N. Equation (^5|) shows that in the limit N — ► oo, 
p is proportional to tit for r s > (37r/4), and equal to N 
otherwise. As expected, the number p of equally occu- 
pied components drops as tit is reduced. Equation (|30| ) is 
equivalent to saying the dimcnsionlcss interparticle sep- 
aration per component for which p = p^ is 




r s (Po 



when Poo/N < 1. The equilibrium value of the N 
MFA energy corresponding to Eq. (|3C|) is 



4\ 2 e 2 



(31) 



3tt J Aireao 



tit- 



(32) 



For arbitrary finite N, Eq. ( p9|) can be used to find 
the interparticle spacing per component ri°\p,p + 1) at 



which the p- and the (p - 
same energies at d = 0: 



r<; 0) (p,p + 1) 



l)-component phases have the 



(33) 




This is the interparticle spacing per component for the 
MFA transition between the p- and (p + l)-component 
phases. It is interesting to compare this result to Eq. j3l| ) 
and to note that 



r(°) (p, p + 1) < r s ( Poc ) < r?> (p-l,p) 



(0), 



(34) 



so that even for finite N, r s (p OCJ ) always gives a value for 
the interparticle spacing that is in the p-component MFA 
phase at d = 0. 

For systems of physical interest containing N/2 layers 
of spin-1/2 particles, it is convenient to work with the in- 
terparticle spacing per layer (rather than the spacing per 
component). This is accomplished by dividing Eq. (p3) 
by \pl. For double-layer systems of spin-1/2 particles 
(N = 4), 



,(o) 



(°)(2,3) = ^(l 



1) w 2.844 
/2/3) w 2.140 



(35) 



(0) 



3tt 



(3,4) = — (^2/3 + 



1.795, 



where the superscript (0) denotes zero layer separation. 
Note that the direct four- to two-component MFA tran- 
sition predicted to occur at 

r s (2, 4) = — (1 + VT72) w 2.011 (36) 

does not exist. In fact, we find that when r s = r s (2,4), 
the gate-balanced (pp — pb) double-layer system is 
always in the three-component MFA phase, regardless 
of the layer separation. The three-component system, 
which at d = has n a -\ — n a [ = = nr/3, has a 
spin-unpolarized subband (n a -\ = rial) with greater den- 
sity than the other subband (n a > rib), which is spin- 
polarized. These features of the three-component phase 
persist at finite layer separations, although the layer im- 
balance is greatly reduced. 



2. Finite layer separation 

Classically (when only the electrostatic energies are 
considered), balanced layers (ni = 712) are obtained when 
Pf = Pb-, in order to make the electric field E12 between 
the inner layers vanish. This result gives the asymptot- 
ically correct behavior for high layer densities and large 
layer separations. At sufficiently low densities and layer 
separations, the exchange energy can dominate the ki- 
netic and electrostatic energies, so that the possibility 
of strong intralayer exchange leading to an interlayer 
charge-transfer instability must be considered. However, 
within the MFA, it can be proved that the inner layer 
densities are always equal, except in the three-component 
phase. 

If the subband densities are equal (n a = rib, the case 
of "pseudospin paramagnetism"), then Eq. (B2) shows 
that ni — ri2- Thus the four-component {n as — jt-t/4) 
and two-component (with n a -\ = nb] = tit/ 2) phases 
have balanced layers. This is because the MFA state 
constructed by occupying equally the single-particle sub- 
band states a and b is equivalent (up to a global pseu- 
dospin rotation) to the MFA state constructed by occu- 
pying equally the single-particle layer states 1 and 2. The 
fact that the ground-state energy in Eq. ( p3| ) is indepen- 
dent of 9 and 4> when n a — nb is due to the invariance 
of the ground-state energy under global pseudospin rota- 
tion. 

If the subband densities are not equal (n a > nb), ex- 
tremizing the total energy per unit area in Eq. ( p3| ) with 
respect to 9 for pf = Pb and t = gives the condition 
sin(26>) = 0. The requirement that the extremum be a 
minimum (i.e., that the second derivative of the total 
energy per unit area with respect to 9 be positive) gives 



if T < 1 
if T > 1 



(37) 



where the interlayer exchange parameter Y is defined in 
Eq. (|2^). The properties of T are described in Sec. A 4 
of the Appendix. Using the inequality 
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e 2 d/2e>V 11 (q)-V 12 (q), 



(38) 



which is true for d > 0, it follows from Eq. (|27j) that for 
d > 0, 



r < 



Lt.L 



1-aas (q) + hbs{q) - 2I abs (q) 



(n a ~ n b ) 2 



(39) 



E 



(n as - n bs ) 2 



(n a - n b ) 



2 ' 



Thus the condition 



{n a \ ~ n b] )(n ai - n H ) > 



(40) 



is sufficient to guarantee that T < 1 for d > 0, so that 
9 = it/2, which balances the layers. The one-component 
phase (n | = ut) satisfies Eq. ([f"o|), so 9 = ir/2 and the 
layers are balanced. No interlayer charge transfer is ob- 
tained in the one-component phase because the combined 
effects of electrostatics and interlayer exchange, which fa- 
vor balanced layers, always dominate the unbalancing in- 
fluence of intralayer exchange: e 2 d/2e + Vi2(q) > Vu(q). 

Within the MFA, the spins in subbands a and b are ei- 
ther completely unpolarized (at higher densities) or fully 
polarized (at sufficiently low densities). Therefore the 
only possible MFA configurations of spin and pseudospin 
which could have unbalanced layers [n a > n b , 9 ^ n/2) 
when pf = Pb would be three-component states with 



M = n ai < ribf, n bl 



0. 



(41) 



We find numerically that such states always have T > 1 
when t = 0. so that sin# = 0. Hence, there is no in- 
terlayer phase coherence and the pseudospins are Ising- 
like. States with n a > n b and sin# = have partially 
unbalanced layers, with the lower-density layer being 
spin-polarized and the higher-density layer being spin- 
unpolarized, even for balanced gates (pp = Pb)- This is 
the behavior we find for the three-component MFA phase. 
For infinitesimal d, one (spin-unpolarized) layer has twice 
the density of the other (spin-polarized) layer, and the 



phase exists in the range r { s (3,4) < r s < r s (2,3). At 
finite d, the three-component phase has only a slight layer 
imbalance and exists only in a narrow region of average 
interparticle spacing around r s r s (2,4). 

The equality of the inner layer densities in the bal- 
anced case (except for the three-component phase) has 
also been shown to be true for the one-component phase 
when intralayer and interlayer oorrelations are included 
within the STLS approximation!^ We note that if inter- 
layer exchange were omitted from the total energy per 
unit area by setting V\2{q) = in Eq. (^3|), then four 
MFA ground states would still be obtained, and the four- 
, three-, and two-component states would be unchanged. 
However, Eq. (^8|) would not be satisfied at small q, and 
at large interparticle distances (at small values of kp, 
or low densities) the vanishing of the second derivative 
with respect to 9 would give the condition cos(2#) > 0, 



so that cos(6*) = ±1. This would produce the inter- 
layer charge instability of Ruden and Wu for the one- 
component phase.lia The fact that the one-component 
phase has equal densities is due to the effects of inter- 
layer exchange. 

Before obtaining the MFA phase diagram for double- 
layer systems, it is convenient to express lengths and en- 
ergies as dimensionless quantities. We therefore express 
the layer separation d and the average density per layer 
np/2 in terms of the effective Bohr radius of the sample, 
ao = AneTi 2 /m*e 2 , where e is the dielectric constant in 
the material, and m* is the effective mass. For GaAs, 
the dielectric constant is e « 13eo- For n-type GaAs, 
m* ss 0.07m e so that ao ~ 9.8 nm, while for p-type 
GaAs, m* « 0.3m e so that ao « 2.3 nm. The average 
density per layer can be expressed in terms of dimension- 
less ratio r s = ro/ao of the interparticle spacing ro for 
a single layer of averaged density ut/2 (defined through 
irr 2 l TiT/2 = 1) to the effective Bohr radius ao- Thus, for 
the same total density, p-type GaAs will have a value 
for r s that is about 9.8/2.3 ~ 4.3 times larger than for 
rt-type GaAs. We define the Fermi wave vector kp in 
terms of the total density tit as 



kp = V 47T71T jp = 



2 2 



pr s a Q 



(42) 



for p = 1,2,4. For the spin-polarized p — 2a state and 
and for the completely unpolirizcd p = 4, kp is equal to 
the Fermi wave vector corresponding to the average den- 
sity per layer (np/ '2). The energy scale associated with 
the effective Bohr radius is vq — e 2 /4neao = h 2 /m*aQ, 
which gives vq m 11 meV for n-type GaAs, and vq « 48 
meV for p-type GaAs. 

The HFA phase diagram for balanced double-layers 
without j-tunncling was obtained by Zheng and co- 
workers,E3 except for the three-component phase. Like 
them, we find that within the MFA, three of the stable 
phases have equal average inner layer densities. Only the 
three-component phase has unequal layer densities. To 
understand the origin of the MFA phases, it useful to 
consider the five terms that contribute to the total en- 
ergy per unit area in Eq. (^J) • The first term of Eq. ( |2^ ) 
is the kinetic (Fermi) energy, which favors distributing 
the particles equally among the subbands and spins. At 
the highest densities, the kinetic energy term dominates, 
and the double- layer system is a four-component spin and 
pseudospin paramagnet: n a j = n a[ = n b ^ = n b i = rip/ A. 
The second term of Eq. (E3T) is the tunneling energy, 
which we take to be zero for now. In general, it favors 
n a > n b (pseudospin polarization), without regard to 
the real spin. The third term of Eq. ( p3| ) is the electro- 
static energy, which vanishes when the gates and inner 
layers are balanced. In general, the electrostatic term fa- 
vors complete screening, which would make n\ = pp and 
^2 — Pb, without regard to the real spin. 

The fourth term in Eq. ( ^3| ) is an intrasubband ex- 
change term that dominates at the smallest densities and 







layer separations. It has the opposite effect of the ki- 
netic energy, eventually producing a one-component spin 
and pseudospin ferromagnet at very low densities and 
small layer separations. The last term (containing the 
interlayer exchange) favors pseudospin paramagnetism 
{n a s — nb s ), but is indifferent to the polarization of the 
real spin, so long as it is the same in both subbands. 
Thus when pp = pb and d > 0, the last term and the 
electrostatic terms are responsible for producing a two- 
component phase that is ferromagnetic in real spin rather 
than in pseudospin: n a ^ = = Uy/2. (At d = 0, two- 
component states that are ferromagnetic in either the 
spin or pseudospin are degenerate. For d > 0, the pseu- 
dospin ferromagnetic p = 2b state is favored only for sub- 
stantial tunneling t and/or layer imbalance \pp —Pb | •) In 
the absence of the kinetic energy term (e.g., in the limit 
to* — > oo), the real spin is always polarized due to in- 
tralayer exchange. However, because of the interlayer 
exchange, the two-component state is always obtained at 
high densities or large layer separations, while the one- 
component state is favored only for low densities and 
small layer separations. 

Finite layer separation (d > 0) differentiates be- 
tween spin and pseudospin, so that the symmetry of 
the problem becomes SU(2) x U(l) rather than the 
CP (3) symmetry at d — 0. At finite layer separation, 
^12(9) < Vi\(q)i an( i the last term in the energy per unit 
area in Eq. (|23j) is minimized by equal subband densities, 
rather than by equal real spin densities. This does not 
effect the p = 1 (fully spin- and pseudospin-polarized) 
or p = 4 (completely spin and pseudospin unpolarized) 
phases. However, when t — and pp = ps, the two- 
component [p — 2) MFA phase has its real spin fully 
polarized and its pseudospin completely unpolarized. 

Finite layer separation also changes the densities at 
which the transitions between neighboring p-component 
phases occur. The densities at which the transitions be- 
tween the phases occur in the MFA can be determined 
by comparing MFA energies. This task is made easier 
by the fact that the MFA always makes the real spins 
in a given subband either completely unpolarized (para- 
magnetic) of fully polarized (maximally ferromagnetic). 
This would probably not be true if correlation-energy ef- 
fects were properly included]!] and it is likely that in 
real double-layer systems, states with partial polariza- 
tion may be stable in some regions of density. The same 
MFA behavior (the restriction to the two extremes of ei- 
ther zero or full polarization) is also found for the pseu- 
dospin when the layers are balanced and the interlayer 
tunneling is zero, except for the three-component phase, 
which has partial pseudospin polarization. 

The four-component and two-component phases both 
have equal subband densities, so their MFA energies are 
independent of the layer separation d. If there were 
a direct transition between these two phases, it would 
be simply a spin paramagnetic to ferromagnetic tran- 
sition in each subband or layer. Therefore, within the 
MFA, such a four-component to two-component transi- 



tion would occur at the same layer density as the spin- 
polarization transition for a single-layer system with a 
layer density equal to the subband densities: i.e., for 
r s (2,4) = 3tt(1 + y/T/2)/8 » 2.011, independent of the 
layer separation. As we discuss below, the direct four- 
to two-component transition is interrupted by a three- 
component phase, which has one subband (layer) spin- 
polarized and the other spin-unpolarized. So it is still 
true that MFA spin polarization transitions occur near 
r s w r s (2, 4). However, the actual value of r s needed 
to obtain spin-polarization in a real sample is likely to 
be significantly higher. For single-layer systems, diffu- 
sion Monte Carlo simulations show that the low-density 
ferromagnetic state predicted by the HFA does occur; 
however, correlation-energy effects move the transition to, 
densities that are probably 100 times lower, to r s w 20x3 
Such high values of r s have been achieved in low-density 
p-type GaAs samples, which possess a larger effective, 
mass (and therefore larger r s ) than n-type samples.EEl 
Large values of the effective mass will favor the existence 
of the lower-component (p < 4) described here, in that 
they increase r s . 

We find empirically that the three-component MFA 
phase has sin 6* = 0. In order to analyze this phase, con- 
sider a MFA ground state with 

n hT = (n T /3)(l + x/2), n H = 0, (43) 
n a ] = n a i = (n T /3)(l - x/4), 

and sin 6* = 0, for < x < 1. Note that the layer imbal- 
ance is given by 

An = n a -n b = (n T /3)(l-x). (44) 

When d = 0, the three-component phase distributes the 
densities equally between the three components (but not 
between the two layers), so that x = 0. As d — > 00, the 
Coulombic cost of layer imbalance becomes prohibitive, 
and x — > 1 so that An — > 0. We plot the layer imbalance 
ratio An/(n T /3) = (1 - x) for r s = r s (2,4) « 2.011 in 

Fig. |. 
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FIG. 2. Layer imbalance ratio An/(riT /3) = 1 — x in the 
three-component phase for r s = r s (2, 4) w 2.011. The imbal- 
ance decreases with increasing layer separation. 



As d -► 0, 



(4/3) (d/a„) 




(45) 



to linear order in d/ao, where the last expression in the 
equation above is for r s = r s (2, 4), which is the only 
value of r s for which the double-layer system is in the 
three-component phase for arbitrary layer separation. As 

d — > oo, 



An 

n-r/3 



1 



oo 3 
d 4 



(1 + l/V2>/2, 



16 d ' 



(46) 



to linear order in a^/d, where the last expression in 
the equation above is for r s = r s (2, 4). Equation fl46| ) 
says that the layer imbalance An in the three-component 
phase is inversely proportional to d/ao for large values of 
d/a . 

For infinitesimal d, the energy per unit area 
of three-component phase increases by the amount 
(e 2 d/8e)(nT /3) 2 to linear order in d, while that of the 
four- and two-component phases are unchanged. Equat- 
ing the three-component energy to the four- and to the 
two-component energies gives 



r.(2,3) 
,r.(3,4) 
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to first order in d/ao as d/ao 
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r s (3,4)«r(°)(2,4) (l- 
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(48) 



to first order in ao/d as d/ao - > oo. Note that both 
r s (2, 3) and r s (3, 4) approach r s (2, 4) in the limit d — ► oo. 
This is because as d — ► oo, the double-layer system con- 
sists of two independent layers, and r s (2, 4) is the in- 
tcrparticle spacing at which the spin-polarized and spin- 
unpolarized energies are equal in a single-layer system. 
Thus, as d — > oo, the energies of the four-, three-, and 
two-component phases are all equal at r s = r s (2, 4). 

We now consider the two-component to one-component 
transition in the MFA. For d = 0, the MFA transition 
to pseudospin fcrromagnetism is equivalent to a real-spin 
paramagnetic to ferromagnetic transition in a single layer 
having total density tit (rather than n-r/2). Thus for 
d = 0, the MFA critical density per layer for the two- 
to one-component transition is exactly half the critical 
density for the four- to two-component transition, so that 
ri 0) (l,2) = 3tt(\/2 + l)/8 w 2.844. By equating the one- 
and two-component phase energies per area, the critical 
density for the one- to two-component transition may be 
obtained: 



ri 0) (l,2) 
r.(l,2) 



where 



z = 2kpd 
= (2 - V2) 



l-(l + l/^)^ ri (.) 
1 - (1 + 1/V2) [1 - S(z)] 

2dV47rn T = AV2d/(r s a ) 
32 d ri 0) (l,2) 



(49) 



(50) 



37r a,Q 

Ti = r(n aT = n T ), and 



S(z) 



dxe~ zx 
2 



3 

2 jo 
3tt 

iz l z 

' 1 - (3tt/32)z+ (l/5)z 2 , z -> 



1 



arccos(a;) — xy 1 — x 2 
M*) - Li(z)] 



(51) 



3tt/4z, 



oo 



Here I\ and L\ are modified Bessel and modified Struvc 
functions of the first kind, respectively. The derivation 
of the above formula is discussed in Sec. A4c of the 



Appendix. Equations ( [49|) and (|50j ) determine r s (l,2), 
which we have plotted as the upper solid line in Fig. |3| 
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d/a o 

FIG. 3. MFA (solid lines) phase diagram for a 
gate-balanced (j>f = Pb) double- layer system with t = 0. 
At the highest densities (smallest r s , bottom of figure), the 
four-component state, which is both spin and pseudospin un- 
polarized, is energetically favored. In a narrow region around 
r s (2,4) « 2.011 (dotted line), the three-component state has 
the lowest energy. At lower densitities, a two-component 
state which is spin polarized but pseudospin unpolarized is 
favored. At the lowest densities (large r 3 , top of figure) a 
one-component state which is both spin and pseudospin po- 
larized is favored. The dotted line at the top is a GRPA 
estimate for the onset of the one-component phase. 

It follows from Eq. @) that r s (l, 2)/ri° } (l, 2) -> 
(1 + d/a ) in the limit d/a — > 0. The value of r s (l,2) 
at very large layer separations (d/ao — > oo) can be 
obtained by setting r s (l, 2)/ri°\l, 2) = in the left- 
hand side of Eq. ([l9j) and then solving the resulting 
equation numerically for z. The result is z — > 4.015, 
or r s (l,2)/ri 0) (l,2) -> 0.495d/a , so that r s (l,2) -> 
lA09d/aQ. We note that if interlayer exchange were 
ignored and the pseudospins treated as Ising variables 
(sin 9 = 0), then the two- to one-component MFA tran- 
sition would occur at ri 0) (l,2)(l + d/ao) (for all values 
of d/ao) and would put all (spin-polarized) particles in a 
single layer. Interlayer exchange makes the layer densi- 
ties equal in the one-component MFA state, and causes 
the two- to one-component MFA transition to occur at 
a somewhat lower value of r s than would be predicted 
using Ising pseudospins. 

B. Infinitesimal tunneling 

In this subsection, we discuss the effects of very small 
interlayer tunneling. The contribution of interlayer tun- 
neling to the MFA energy per unit area is given by 
—2t(n a — rib) sin 9. Finite interlayer tunneling thus has 
two important effects. First, it removes the Ising charac- 



ter of the spins by making sin 9 > 0. We discuss this 
in more detail in the next subsection. Second, it al- 
ways produces some degree of pseudospin polarization 
(n a — rib > 0). We have parametrized the dependence of 
the subband splitting (n a — rib) on t in terms of the pseu- 
dospin Stoner interaction parameter, which we calculate 
in the second subsection below. 

1. Effect of tunneling on phase angle 

We first consider the effect of interlayer tunneling on 
the interlayer phase angle 9. Extremizing the energy per 
unit area [Eq. (23)] with respect to 9 for t > gives two 
possible solutions. The solution 



cos 9 = 



(52) 



is always an extremum, and is a minimum whenever T < 
1, which is true for all phases except p = 3 near balance 
(pf ~ Pb)- Of course, if t were large enough to cause full 
pseudospin polarization (n a — tit), then cos 9 = 0. 
If r > 1 (three-component phase) then the solution 



t 



suit 



(r - l)e 2 d(n a - n b )/2e 



(53) 



gives the correct minimum, provided that the right-hand 
side is positive and does not exceed unity. Even in the 
three-component phase, t will increase (n a — rib), but 
since (n a — rib) is already finite at t = 0, sin# will be 
proportional to t, to first order in t. 

The effect of including a small amount of interlayer 
tunneling can often be described perturbatively, and pro- 
duces a smooth increase in the pseudospin polarization 
(n a — rib) that is proportional to t. However, the re- 
sults obtained in the combined limits t — > and d — > 
can depend on the order of these limits. For exam- 
ple, at d = 0, the lowest-energy two-component phase 
has n a i = n a i = n r/2 (pseudospin-polarized but spin- 
unpolarized) for arbitrarily small but finite t. How- 
ever, at t = 0, the lowest-energy two-component phase 
has n a -\ ~ n fcT = nr/2 (spin-polarized but pseudospin- 
unpolarized) for arbitrarily small but finite d. By com- 
paring the MFA energies per unit area of two competing 
p = 2 ground states with cos 9 = (spin-polarized p = 2a 
versus pseudospin-polarized p = 2b), it can be shown 
that the pseudospin-polarized two-component (p = 2b) 
ground state requires 



t > 



e 2 driT T ^, . . e 2 driT 



-r(n aT = n ai = n T /2) = 



r x (p = 2) (54) 



e 2 driT 16 



[l-S(z)} 



3irz 
e 2 driT /16e, 



k^d 



(4/3V27r)e 2 v /^/47re, k F d 





oo. 



Here 



2k F d = 2d^/2nriT, and S(z) is defined in 



Eq. @. Therefore it is the size of t relative to e 2 driT/16e 
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that must be considered as both t and d approach zero 
in the two-component phase. A rough estimate of the 
minimum tunneling energy t c necessary to obtain the 
pseudospin-polarized two-component phase may be made 
by calculating e 2 driT/16e for the smallest value of rip 
that still gives the two-component phase: i.e., for r s ~ 
A 0) (l,2) 2.844. For d = 10 nm, this gives t c ~ 0.7 
meV for n-type GaAs and t c ~ 13 meV for p-type GaAs. 
The differences between these two values of t c arise from 
the fact that n-type GaAs spin-polarizes at a much lower 
density than p-GaAs, due to the differences in the effec- 
tive mass (and therefore in r s .) We stress that there ex- 
ist two p = 2 two-component states, which can be either 
spin- or pseudospin-polarized. The p — 2 state that we 
focus on most will be the spin-polarized state (p = 2a) . 

2. Pseudsopin Stoner parameter 

The difference An = n a — rib between the subband 
densities (obtained from SdH measurements) of a bal- 
anced double-layer system is often used as a measure 
of the size of the interlayer tuncling matrix element t 
by applying the formula n a — ft*-j= (p/2)tvQ, which is 
valid for noninteracting particles.cJ Here p = 4 for spin- 
unpolarized particles, p = 2 for spin-polarized particles, 
and Vq = m* / (n/h 2 ) is the two-dimensional density of 
states per unit area. For interacting electrons or holes, 
the pseudpspin Stoner interaction parameter / is defined 
througbEa 



1-i 



7T/2 A _ -2fc F dsin0\ 
d0 ~ 



n a - n b 



(p/2)tv 



1 



(55) 



For noninteracting particles, 1 = 0. For interacting par- 
ticles, / is a function of r s and d/an; in general, it also 
a function of t. We do not consider p = 3 or p = 1, 
since both the three- and the one-component phases have 
n a > rib at t — 0, corresponding to I = 1. The onset of 
SILC (n a > rib and sin 9 > even when t = 0) occurs at 
the two- to one-component transition, and corresponds 
to / — > 1 for p = 2. 

The Stoner interaction parameter I can be calculated 
analytically in the limit of vanis hing interlayer tunnel- 
ing (t — ► 0), as we show in Sec. |A 5| of the Appendix. 
(For finite t, we calculate I numerically.) The basic idea 
is to start with equal subband densities (n n = rib) and 
infinitesimally small t, and then calculate the change in 
energy due to moving an infinitesimal amount of charge 
from subband b to subband a. Minimizing the change in 
energy with respect to the amount of charge transferred 
between the subbands gives the linear response to small 
interlayer tunneling and yields the following expression 
for the t — > limit of the Stoner interaction parameter: 



— | 1 — p— kpdTn 

2irek F V 1 4 ° 



2Vn(2fc F ) 



dx 



(56) 

[Vn(2k F x) - V 12 (2k F x)} 



where kp = y/ Awrvp /p for p = 2,4, and r n = r(n a — > 
rib). Equation ( |56| ) is equivalent to the generalized 
random-phase approximation (GRPA) result in Eq. (14) 
of Ref. ||. 

Note that 1=1 (when p = 1) corresponds to SILC, 
and that the GRPA result for the phase boundary, which 
is shown as the dashed sloped line in Fig. ||, is different 
(has larger r s ) than the MFA result. That the GRPA 
gives a higher value for r s (l,2) is not surprising, given 
that the GRPA goes beyond MFA and contains correla- 
tion effects in an approximate way. To lowest order in 
d, 1(d)/ 1(0) = (1 - d/ao), so that for d -> 0, the GRPA 
gives a higher critical value of the interparticle spacing 
for SILC than the MFA: r s (l,2) w tt(1 + d/a Q ). As ex- 
pected, a similar calculation of the linear response of the 
real spins to a weak Zeeman field shows that a hypotheti- 
cal four- to two-component GRPA transition would occur 
at twice the density of the d = two- to one-component 
transition, i.e., at r s (l, 2) = Tr/y/2. However, as with the 
MFA, we expect that a three-component GRPA phase 
preempts any direct four- to two-component GRPA tran- 
sition, and that a three-component is always obtained 
within the GRPA at r s = ttA/2. 

As pointed out in Ref. E3, interactions enhance the 
subband splitting (/ > 0) for kpd < 1.13, but reduce 
the splitting (/ < 0) for kpd > 1.13. The critical 
value for r s at large d that separates the one- and two- 
component phases is determined by solving Eq. (pq) for 
1 = 1 and r s — > oo. Asymptotically (for r s ,d/ao — > 
oo), it occurs at the same value of r s that has 1 = 

i.e., k F d -> 1.134 or r a (l,2)/r£ 0) (l,2) -> 0.620d/a , 
which gives r s (l,2) -> 1.764d/a . Within the GRPA, 
SILC occurs once the interparticle spacing is roughly 
twice the interlayer spacing. In the limit kpd — > oo, 
1(d)/ 1(0) = 1 - (l/2)[ln(4A;FcO +7], where 7 w 0.5772 
is Euler's constant. In Ref. |2^ it is argued that although 
1(d) is large and negative as d — > 00, its apparently di- 
vergent behavior is an unphysical artifact of the GRPA. 

We have calculated the Stoner interaction parameter 
I for a few values of the interlayer tunneling t in Figs. ^ 
(n-type GaAs) and || (p-type GaAs), for a hypothetical 
sample with total density nx = 10 n cm~ 2 . The com- 
plete polarization of the pseudospin (n a = tit) is in- 
dicated by the mesa-like regions where / becomes flat: 
n a ^ — n a i — n T/2 for n-type GaAs and n a? = ny 
for p-type GaAs. Increasing the size of t favors pseu- 
dospin polarization and allows it to persist to larger val- 
ues of kpd. Note that at fixed t, the Stoner interac- 
tion is equal to I max = 1 — 2<^n/ n T as kpd — > 0, so 
that / max decreases with increasing t. The fact that 
Ima,x(t = 0.1) > I max (t = 0) is an artifact of the order in 
which the limits t — > and d — > are taken: for any finite 
t, the two-component phase will have n a j = n a i = rip/2 
as d — -> 0, but for any finite d, the two-component 
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phase will have n a f = n&f = n-r/2 as f — > 0. For 
n-type GaAs with t = 1 meV (dotted line in Fig. 
the system goes through three phases as a function of 
kpd: (1) a pseudospin-polarized two-component phase 
(fiat = n <4 = n T/2) for fc^d — ► 0, (2) a three-component 
phase (n a | = n a | < ntf) for intermediate values of fc^d, 
indicated in by the "missing piece" on the right side of the 
mesa in Fig. |f|, (3) a real-spin-polarized two-component 
phase (n a | = = np/2) for larger kpd — > oo. 



1-0 pr 
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FIG. 4. Stoner interaction parame- 

ter 1=1 — 2tvo/(n a — rib) versus kpd for n-type GaAs at 
total density n T = 10 n cm" 2 for t = 0,0.1, 1.0 meV. 

SILC is indicated by having I(t = 0) = 1 and occurs 
only for p-type GaAs (Fig. |, for kpd < 0.7): I(t = 0) is 
directly proportional to vo, and therefore to the effective 
mass of the particles, so that SILC is more likely to be 
observed in p-type (to* /m e 0.3) rather than n-type 
(m*/m e ps 0.07) GaAs. 
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FIG. 5. Stoner interaction parame- 

ter 1 = 1 — 2tV(j/{n a — n ) versus kpd for p-type GaAs at 
total density n T = 10 n cm" 2 for t = 0, 0.1, 1.0 meV. 



V. EFFECT OF BIAS 

In this section, we study the effect of bias (pp ^ pb, 
due to applied gate voltages) on the subband and layer 
densities. The classical results 

ni=p F <d(pp), u 2 =Pb&{pb), (57) 

and 

V F = eDpEp + edE 12 + eV F 0) (58) 

give the asymptotically correct behavior for high layer 
densities and large layer separations. Double-layer sys- 
tems at low densities and small layer separations show 
measureable deviations from the classical behavior, most 
notably because of quantum-mechanical exchange. We 
shall find it convenient to study the effects of layer im- 
balance by fixing the total density (np = pp + pb) and 
then varying the gate-imbalance parameter £, defined by 

C = — VB = {pf -PB)/n T . (59) 
Pf +Pb 

The case of balanced gates (pp = pb) corresponds to 

C = o. 

In the presence of bias and/or tunneling, there are 
five possible noncrystalline MFA ground states, which we 
write in order of increasing r s in Table | Only the last 
(p = 1) phase can exhibit SILC. The last two phases are 
pseudospin-polarized and can have < sin# < 1, with 
the value of 9 being determined by the layer separation, 
density, bias, and tunneling. The first three phases are 
Ising-like (sin# = 0) when t — 0, and therefore do not 
involve interlayer exchange, in the absence of interlayer 
tunneling. Note that there are twop = 2 phases, one that 
is pseudospin-polarized but spin-unpolarized (p = 2b), 
and another that is spin-polarized but pseudospin un- 
polarized (p = 2a). The pseudospin-polarized p = 2b 
requires bias and/or tunneling. For all MFA phases, the 
real spin is either fully polarized or completely unpolar- 
ized. 



TABLE I. Subband occupancies of the five possible 
non-crystalline MFA ground states for a double layer system 
with gate imbalance (\pp — Pb\ > 0) or tunneling (i > 0). 



p = 4 


rial = n <U ^ n f>T = n H > 




p = 3 


n a -\ = rial >ribi>Q 


nbi = 


p = 2a 


rial >n b i>0 


rial = nb\. = 


p = 2b 


rial = n al > 


n b i = ribi = 


p=l 


rial = n T 


Hal = n bl = n bl = 
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We begin this section by studying the case of van- 
ishing interlayer exchange, which is the relevant situa- 
tion for the majority of double-layer samples that have 
been studied experimentally, except in the quantum Hall 
regime. tl This is the case when layer separations are suffi- 
ciently far apart that interlayer exchange is negligible in 
zero or weak magnetic fields. A model without interlayer 
exchange is able to fit existing double-layer data well, 
and obtains the correct four-, three-, and two-component 
phases, although it fails to properly describe the one- 
component phase. We then include the effects of in- 
terlayer exchange and solve for the subband and layer 
densities using our MFA, which allows for the possibility 
of interlayer exchange, even in the absence of interlayer 
tunneling. Finally, we give a full treatment of the one- 
component phase within the MFA, including bias and 
tunneling. 

A. No interlayer exchange 

It is simplest to begin our study of the effects of layer 
imbalance {j>f ^ Pb) by first considering the limit of van- 
ishing interlayer exchange. This limit is relevant to most 
of double-layer samples that have been studigd experi- 
mentally, except in the quantum Hall regime.cl It corre- 
sponds to layer separations that are sufficiently far apart 
that interlayer exchange is negligible in zero or weak mag- 
netic fields. We shall also demonstrate that a very sim- 
ple model which assumes that the particles are always 
spin-unpolarizcd gives a good fit to existing data on the 
subband occupancies of double-layer systems, except at 
low densities. 

1. No tunneling 

Interlayer exchange is negligible when the layer sep- 
aration is large compared to the interparticle spacing. 
This condition may be expressed in various ways, e.g., 
kpd > 1 or r s < d/a , and is satisfied for most samples. 
In this limit, we ignore interlayer tunneling and correla- 
tions and write the exchange-correlation energy in Eq. (||) 
as 



e(m, n 2 ) « e(m) + e(n 2 ), 



(60) 



where s(n) is the sum of the kinetic (Fermi), exchange, 
and correlation (but not the electrostatic) energies, for a 
single-layer two-dimensional electron gas (2DEG) of den- 
sity n. In the absence of interlayer exchange and tunnel- 
ing, we can work directly with the layer densities rather 
than the subband densities because the pseudospins are 
Ising variables: 



n a = max(ni,n 2 ), n& = min(ni, n 2 ), 
= 7T<d(p B -Pf) 



(61) 



7T, p F < p B 
0, p F > Pb 



No interlayer exchange is required to correctly describe 
the four-, three-, and two-component MFA phases at t = 
0, since all have sin (9 = 0. Our analysis of these phases 
therefore proceeds as before. It is nonetheless useful to 
look at the phases of system in terms of the equilibrium 
and stability conditions of Sec. ||, which we do below. 

Although it would be straightforward to include in- 
tralayer correlation-energy effects in s(n), we do not do 
so here for the sake of simplicity. Even so, the resulting 
approximate model does a very good job of fitting SdH 
data until the layer densities get so low that they violate 
our initial assumption that kpd is large. We therefore 
begin with 



e(n) 



E 



3\/lT 47T6 



3/2 



(62) 



n 2 /2is Q - (8/3^/2^) (e 2 /47re)n 3 / 2 , n > n c 
n 2 /v - (8/(3V^)(e 2 /47re)n 3 / 2 , n < n c ' 

where n c is the critical density for the MFA spin- 
polarization transition in a single-layer system. The con- 
dition n > n c corresponds to spin-unpolarized electrons 
and occurs at higher densities, whereas n < n c corre- 
sponds to spin-polarized electrons and occurs at lower 
densities. When sin^ = (which is the case we are con- 
sidering in this section), n c is the critical density of the 
MFA (or the HFA) spin-polarization transition for a sin- 
gle layer, which occurs when the single-layer r s has the 
value r s (2,4) = (3tt/8)(1 + y/T/2) « 2.011, shown as 
the dashed line in Fig. |[ Note that r s = 1/ y/ una 2 
for a single-layer system of number density n, which 
gives n c al = (2/3)(4/tt) 3 (1 - 2\/2/3) « 0.07870. In 
the MFA, which is equivalent to the HFA for balanced 
layers, the spin polarization is either completely unpo- 
larized (at higher densities) or completely polarized (at 
lower densities). Correlation-energy effects probably pro- 
duce a range of intermediate spin polarizations. 

The chemical potential measured relative to layer i is 
Hi = fi(rii), where 



fi{n) = de(n)/dn 

n/v Q - (4/V27r)(e 2 /47re)V^, n > n c 
2n/u - (4/y/^)(e 2 /4ire)y/n~, n < n c 



(63) 



where we have used Eqs. (||) and (^2|). The values of the 
layer densities can be determined by using Eq. ( |63|) in 
the equilibrium condition of Eq. (Q), fj,\ — [ii = eE^d. 

The electronic lengths Sy that determine the Eisen- 
stein ratio Re [Eq. (pi])] and the condition for stability 
against interlayer charge transfer [Eq. (0)] can be calcu- 
lated from Eqs. (jfy and (^). Ignoring interlayer cor- 
relations as in Eq. (|60|), gives sy = for i ^ j, and 
Si = su = s(ni), where 



s(n) 
a 



e a 



On 



so that sin 6 = 0. 



1/4 - (y/2/ir)/(iira \/n), n>n 

1/2- (2/vW(47raovH n<n 
1/4 - (V2/47r)r s , n > n c 

1/2 - (l/27r)r„ n < n c ' 
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Note that the MFA compressibility k can be calculated 
from the above result using Eq. (||), and that s(n) and 
therefore k are negative at sufficiently low densities. The 
length s(n) and the compressibility k jump discontinu- 
ously at n — n c , where the ground state changes abruptly 
from spin-unpolarized to spin-polarized. For densities 
just above n Cl s(n+)/ao ~ 0.0237, while for densities 
just below n c , s(n~)/ao ~ 0.1799. Sensitive measure- 
ments of the interlayer capacitance (e.g., the Eisenstein 
ratio Re) could detect the exchange-driven spin polariza- 
tion of a 2DEG through its effect on the compressibility, 
especially in p-type GaAs samples when the density in a 
layer could be made small enough to polarize the holes 
in that layer. 

It is straightforward to calculate the effect of spin po- 
larization in a low-density layer on Re using Eqs. (113) 
and (|64|). For the usual case in which the interlayer 
separation d is substantially larger than the electronic 
lengths s(n,i) (i.e., for d/ao 3> 0.2), the MFA gives an 
abrupt jump in Re by almost a factor of 8, from approxi- 
mately 0.0237a /<i for densities just above n c , to approxi- 
mately 0.1799a /<i for densities just below n c . Of course, 
the MFA overestimates the size of the jump, but it is 
nonetheless plausible that measurements of Re could de- 
tect changes in the compressibility of a low-density layer 
due to the exchange-driven polarization of the spins. 

2. Unpolarized spins 

Ruden and Wu assumed not only that the pseudospins 
were Ising-like (sin 9 = QL but also that the real spins 
were always unpolarized.E£l This limited the phases they 
found &tpE = Pb to two: pseudospin- unpolarized (p = 4) 
at high density and pseudospin-polarized (p = 2b) for low 
density. It is straightforward to compare the energy of 
the four-component phase with that of the hypotheti- 
cal pseudospin-polarized (p = 2b) phase of Ruden and 

Wu and show that they are equal when r s /ri°\l,2) 
1 + 2d/ciQ. Although neither assumption was, strictly 
speaking, correct, it is an interesting and useful fact that 
making such assumptions can yield a simple model that 
fits experimental data for layer densities versus gate bias 
quite well, except at thaJowest densities. Figure ^ shows 
experimental SdH dataEj and a theoretical fit from a sim- 
ple theory that ignores interlayer exchange and takes the 
spins to be unpolarized. The value of the interlayer sepa- 
ration d used in the model is taken to be a fitting param- 
eter. The values of d that we obtain with this simplified 
model always locate the idealized two-dimensional elec- 
trons layers inside the confining quantum wells, although 
d always seems to be somewhat larger than the midwell 
to midwell distance. 
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FIG. 6. Plot of the experimental (circles and triangles) 
and calculated (solid curves) layer densities n\ and Ti2 versus 
front-gate number density pf- A simple model which neglects 
interlayer exchange, intralayer correlations, and spin polariza- 
tion and assumes zero layer thickness yields a good fit to the 
experimental densities, except when n\ becomes small. 

Note that the experimental data fits very well almost 
everywhere, except where the density n\ is very small 
when layer 1 is near depletion. Here the simplest the- 
ory (which omits the possibilities of spin and pseudospin 
polarization) erroneously predicts an interlayer charge 
transfer instability. As we discussed in Sec. |J, such 
an instability is unavoidable when interlayer correlations 
(S12) are neglected. According to the stability criterion of 
Eq. (|7|), the interlayer charge-transfer instability should 
occur in the simplest theory (spin-unpolarized, no inter- 
layer exchange) when 

r sl > it V2(l + 2d/a ) - r s2 , (65) 

where r S i = 1 / \j 'imia^ is the dimensionless interparti- 
cle separation in layer i, and we have taken ri\ < 112 (or 
r si > r s2)- At balance, setting r sl = r s2 in Eq. ( p5| ) gives 
r s = (7r/v / 2)(l + 2rf/a ) as the GRPA version of the crit- 
ical particle separation for Ruden and Wu's hypothetical 
pseudospin-polarization (p = 4 — > p — 2b) transition. 

[In the MFA, r s (26,4) = ^ 0) (2,4)(1 + 2d/a ) gives the 
Ruden- Wu hypothetical pseudospin-polarization transi- 
tion at balance.] 

Note that even in the limit of equal layer densities and 
zero layer separation, r s \ > \/tt/2, which is the GRPA 
value of the particle separation for a single layer to spin 
polarize. Thus even in a theory that neglects interlayer 
correlations, the particles in the lower-density layer (or 
both layers, if they are balanced) spin-polarize before the 
layer empties out. As noted in Ref. [l7], the spin polariza- 
tion of the electron*, predicted by the HFA was ignored 
by Ruden and Wu.Ej However, including the spin polar- 
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ization does not eliminate the interlayer char ge t ransfer 
instability, which according to Eqs. (0) and (|64]) would 
occur at 

r sl > 2tt(3/4 + d/o ) - r s2 /\/2, (66) 

when layer 1 is spin-polarized but layer 2 is not, and at 

r sl >2ir(l + d/a )-r s2 , (67) 

when both layers are spin polarized, if interlayer correla- 
tions could be ignored. 

3. LDF model 

We now introduce a tight-binding local density func- 
tional (LDF) model, which includes the effects of inter- 
layer tunneling in a simple way. We shall follow the pre- 
vious section and ignore for now the effects of interlayer 
exchange, and we shall even take the electron densities to 
be always spin-unpolarized. Such an elementary model 
is capable of fitting experimental data quite well, despite 
its simplicity. 

The Kohn-Sham single-particle equations for our tight- 
binding LDF model is conveniently expressed as a 2 x 2 
matrix equation: 



where 




(A) 

H 



where 



(68) 



(69) 



represents the "on-site" energy of layer j, and the tun- 
neling matrix element —t is off-diagonal. The amplitude 
of the wave function for subband A (A = a, b) in layer j 
is , and the subband energy is E\ . The Hartree con- 
tribution to the on-site energy enters via the interlayer 
electric field E\2, as shown in Eq. (Q). The intralayer 
exchange and correlation contributions to the on-site en- 
ergy are given by the exchange-correlation potential fj, xc . 
In LDF theory, [i xc {n) is equal to the derivative, with re- 
spect to density, of the exchange plus correlation energies 
per unit area of a two-dimensional single-layer system of 
uniform areal density n. Equivalently, Hxci^j) is equal 
to /ij [see Eq. (^)] minus the kinetic energy contribution 
to jjij . For simplicity, we do not include intralayer corre- 
lation energy contributions to fi xc , so we write 



4 e 2 

Hxc(n) « =- — spn. 

y/ir 47TC 

The density in layer j is given by 



= 5>a|. 



,W|2 



(70) 



(71) 



N x = (E F - E x )u @(E F - E x ) 



(72) 



is the areal-density contribution from subband A, and 
= m*/irh 2 is the two-dimensional density of states 
for noninteracting particles. The self-consistency of the 
Kohn-Sham equations enters via Eqs. (71) and (|72j) , 
since the layer densities rij, together with the gate den- 
sities p a , determine the interlayer electric field E12 ap- 
pearing in Eq. (|69|). The Fermi energy E F is chosen so 
that the sum of the subband densities N\ is equal to the 
total density p F + Pb- 

This simple LDF model, which takes the layers to have 
zero thickness and assumes that the real spins are unpo- 
larized, is capable of fitting the experimental layer den- 
sity data closely. This is illustrated in Fig . 0, which shows 
SdH data taken from sample A of Ref. EO. The front-gate 
voltages used in Fig. were calculated using Eq. (|5§|). 




A=l 



V F (V) 

FIG. 7. Experimental subband densities for n a (squares) 
and rib (circles) versus front-gate voltage, together with fit 
(solid curves) from tight-binding LDF model. 



B. Interlayer exchange 

We now allow for the possibility of interlayer exchange 
in biased systems. We found that for balanced systems, 
interlayer exchange becomes important only at low den- 
sities and small layer separations, occuring only in the 
one-component phase. In this section, we explore the ef- 
fect of interlayer exchange on biased double layers, and 
find that it can reduce or suppress interlayer (although 
not intersubband) charge transfers. We find that bias 
always increases pseudospin polarization and sometimes 
reduces the total density required to achieve SILC. 
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1. Interlay er phase angle 

When the interlayer tunneling t is zero, it is possible to 
determine the equilibrium value of 8 that minimizes the 
total energy per unit area, in terms of the equilibrium 
subband densities. Extremizing Eq. ( |23| ) with respect to 
8 for n a > rib gives 



cos 8 



X 




(73) 



1 



where T is defined in Eq. (p7|). In the four-component 
phase, the kinetic energy dominates over the exchange 
energy, and (n a — n^) < \pp — pb\, so that sin 6 = 0. In 
practice, the only time that we find | cos 8\ < 1 (for t = 0) 
is in the one-component phase. For Eq. ( J73| ) to minimize 
the energy with respect to 8, the second derivative of 
the energy with respect to 8 must be positive, which is 
equivalent to requiring that Y < 1. If T > 1 (which is 
true for the p = 3 phase near balance), then cos 6* = ±1. 
It follows from Eqs. (§2|) and d73l) that 



(m - n 2 ) 



~(n a - rib), X < —1 

(PF-PB)/(i-r), -i<x<i 

(n a - rib), X > 1 



(74) 



In the special case of interest where \X\ < 1 (which re- 
quires p = 1), we can calculate the Eisenstein ratio Re 
at the point where the layers are balanced: 



Re(pf = Pb) 



1 r 

-^i 1 - SpbISpf)-^—^' 



(75) 



where we have used Eqs. (10) and (|74j). (Recall that 
Sps/SpF ~ if the back-gate voltage Vb is kept fixed, 
but that Spb/Spf — —1 if the total density is kept fixed.) 

Within the tight-binding model of tunneling, the MFA 
model shows that any finite value of the tunneling matrix 
element t prohibits either layer from completely emptying 
out, regardless of the gate charges (j>f,Pb)- Extremizing 
the energy per unit area [Eq. (p3|)] with respect to 8 for 
t > and pjf / pb shows that neither sin 8 = nor 
cos 8 = are ever local extrema. When t is small and 
pp « pb, including any negative values of pp (i.e., for 
£ < — 1, roughly speaking), then extremizing the energy 
per unit area yields 



t/(e 2 dn T /2e) 

[ici-a-nr 



(76) 



to lowest order in t/(e 2 dn T /2e). This result is a local 
minimum provided that it is positive, which is true in the 
limit we are considering here. It follows from Eq. ( [22] ) 
that ni/riT ~ (<fi/2) 2 , so that n\ cx t 2 . Hence n\, al- 
though small, is always nonzero for t > 0. In actual sam- 
ples, large bias changes the tunneling matrix element t. 



Sufficiently large bias shifts the bottom (minimum en- 
ergy) of the wells relative to one another so greatly that 
t can be driven (for all practical purposes) to zero. 

Interlayer exchange is significant only when the layer 
densities and their separation are sufficiently small. In 
order for interlayer exchange to contribute, there must 
be (1) more particles in one subband than another (n a > 
rib), and (2) nonzero sin# (8 ^ 0,ir). Thus, for example, 
the case of balanced layers with 8 — ir/2 at very low 
densities (n a | = rip) is a situation in which interlayer 
exchange contributes strongly. We discuss this case in 
Sec. VI, and we find there that interlayer exchange does 
indeed suppress interlayer charge transfer. For the case of 
unbalanced layers at high total density, generally n a > 
rib] but when t = 0, 8 is usually equal to or 7r, so 
that in the MFA, interlayer exchange does not contribute. 
Near depletion, where one of the layers empties out, the 
situation is not as clear, so we now analyze that situation 
in some detail later below. 

Even with interlayer exchange, it turns out that the 
MFA model is always unstable with respect to an abrupt 
exchange-driven intersubband charge transfer from (low- 
density) subband b to (higher-density) subband a when 
the particle density in subband b gets small enough. The 
abrupt intersubband charge transfer is probably an un- 
physical feature of the MFA model that is not observed in 
real experiments. We believe that a proper treatment of 
the correlation energies (which have been entirely omit- 
ted here) would help fix this shortcoming. Nevertheless, 
we can still investigate what the MFA has to say about 
interlayer charge transfer, since the subbands are in gen- 
eral different from the layers, when we include interlayer 
exchange. 

2. Subbands densities nearly equal 

We now consider the limit in which An = (n a — rib) <§C 
tit, so that the double-layer system is only slightly 
pseudospin-polarized. This will be the case for the two- 
(p = 2a) and four- (p — 4) component ground states for 
small t and \pf~ Pb\- We begin by computing the change 
in the ground-state energy per unit area due to chang- 
ing An from zero to a small but finite value. Expanding 
Eq. (p3|) to second order in An gives 



A£ _ (An) 2 



L X Ly 

e 2 d 



pvo 



t An sin 8 — 



(An) 



47T6 ^JpTTTlT 

e 2 d ( An 



(77) 



+ — [An cos 8-{p F - pb)] + T 



where 



To 



lim F. 

ia—>nb 



(78) 



The quantity Tq is calculated in Sec. A4b of the Ap- 



pendix. Extremizing Eq. ( [77] ) with respect to 8 gives 
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-2t cos 9 



~2e 



(pf - Pb) sin# 



(79) 



(1 



Tn) An sin 9 cos ( 

; 2e 



By taking the second derivative of Eq. (77) with respect 
to 9, we find that Eq. ( f79|) is a local minimum provided 
that 



2tsin9 



-(1 

Except for p - 



e a 
.e 2 d 



Pb) cos 9 



(80) 



2c 



Ancos(26>) > 0. 



3 near balance, < To < 1, so that 

o<i-r <i. (8i) 



Hence Eqs. (|79|) and (gOJ) imply that if t > but \p F - 
Pb\ = 0, then cos 9 = so that 9 = ir/2 (except for 
p = 3.) If, however, \p F — p B \ > but t = 0, then 
s'm9 = so that 9 = 0,ir (except for p = 1.) There is 
thus a competition between the effects oft and \pf~Pb\- 
If neither t nor \pp — pg\ is zero, then sin# ^ and 
cos(? ^ 0. In the limit that (e 2 d/2c)\pF — Pb\ <C t, then 



{e 2 d/2c){p F -p B ) 
2t 



< 1 



(82) 



for An — > 0. On the other hand, in the limit that t <C 
(e 2 d/2e)|p F then 



2t 



(e 2 d/2e)[\p F -pB\-(l-T )An} 



< 1. (83) 



In general, 9 must be solved numerically. 

Extremizing Eq. uTA) with respect to An gives 



An = 



2t sin# 4- (e 2 d/e)(pF — Pb) cos9 

4 



^(cos 2 6» + r sin 2 1 



(84) 



4-7re y/jmrTr 



which is a local minimum provided that its denomina- 
tor is positive, as may be seen by computing the second 
derivative of Eq. (^) with respect to An. Note that both 
interlayer tunneling (t) and gate bias (\pf — Pb\) produce 
pseudospin polarization (increase An.) 

When pf = Pb and t is (arbitrarily) small but nonzero 
so that 9 = 7r/2, Eq. ( |S4| ) yields the pseudospin Stoner 
enhancement factor in Eq. ([56|) . This is discussed in more 
detail in Sec. |IVB2| . 

When sin 9 = (which requires t — 0), then Eq. j8^ ) 
gives 



A? 



d/ao 



\'Pi- 



Pb\ 



d/a + (2/p)[l-r s /(7T^2/p)} 



(85) 



This is the case for the p = 2a and p — 4 states for 
An/riT -C 1. It follows from this that near balance (pf ~ 
Pb), the Eisenstein ratio for fixed ps is given by 



Rr — — 



(p/2)d/a 



1 - r s /(7r A 



(86) 



where we have use of Eq. (|1C|). Equation (|35|) says that, 
when sin# = 0, then An < \pf — Pb\ for small An/nr, 
provided that 



r s < 7Ta 



(87) 



Now, the GRPA estimate of the interparticle separation 
required for spin polarization in a single two-dimensional 
layer is r s = n/^/2, which is just the right-hand side of 
Eq. ( |S7j ) for p — 4. We therefore expect that, as long 
as the ground state has four components, it will be true 
that An < \pf — Pb\, and this is indeed what our MFA 
calculations find for small r s . 

According to Eq. (|8^), An < \pf~ Pb \ for small An/n^ 
for p = 2a until r s > it. The interparticle separation 
r s = 7r is also the GRPA estimate of (1, 2) required for 
pseudospin polarization. Thus for d > we expect that 
at higher densities (small r s ), An < \pf~Pb \ throughout 
the p = 4 state and in the low-r s region of the p = 2a 
state, but that An > \pp — Pb\ for the high-r s region of 
the p — 2a state, at least for small An/nr- This is in fact 
what our MFA calculations show. It is also true that in 
the p = 3 phase (which has An > even when pp = Pb), 
An > \pf~ Pb \ for C ^ 1> although not for £ on the order 
of one. Of course, for r s sufficiently small, pseudospin 
polarization occurs so that An = rip > \pf — Pb\- 

Figure || shows a plot of the subband densities n as 
and the ratio (pp -~PB/(n a — rib) versus r s for fixed layer 
separation d/ao — 5 and fixed layer imbalance £ = 0.2. 
For small r s , the p = 4 phase is obtained, and (pf — 
PB)/(n a — ub) > 1. For r s fa 2.011, the p — 4 phase is 
obtained, and (j>f — PB)/(n a — tlb) < 1- For larger r s , 
the spin-polarized p — 2a phase is obtained. Note that 
for p = 2a the ratio (pf — PB)/{n a — ns) is larger than 
one for smaller r s , but smaller than one for larger r s . 
For even larger r s (not shown), the p = 1 phase would 
be obtained with (pp — Pb)/ {n a — tib) = C which in this 
case has £ = 0.2. 
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FIG. 8. Normalized subband densities and the ratio 
(pf —pB)/{n a - fig), versus average interparticler separation 
r s , for fixed d/ao = 5 and fixed gate imbalance parameter 
C, — 0.2. The p = 4, p — 3, and p = 2a phases are obtained 
successively as r s is increased. 



If the denominator in Eq. (84) is not positive, then 
the global minimum for the energy per unit area occurs 
for An = tit, corresponding to pseudospin polarization. 
Thus the GRPA condition for stability against abrupt 
intersubband charge transfer is just the condition that 
denominator in Eq. (p3) be positive, or equivalently 



d 
«o 



(2/p)[1 



(cos 2 9 + T skr 



> 0. 



When sin 9 = (which requires t = 0) so that the pseu- 
dospin is Ising-like and interlayer exchange does not con- 
tribute, then Eq. ( pq ) is equivalent to the stability con- 
dition against abrupt interlayer charge transfer given in 
Eq. (0), when the electron lengths s(n) are approximated 
by Eq. (|64|). When cos9 = (e.g., when pp = pb and t 
has any finite positive value) , then the violation of the in- 
equality in Eq. (^8|) is equivalent to the condition that the 
pseudospin Stoncr chanccmcnt factor /, given in Eq. (56), 



is equal to one, which signals the transition to pseudospin 
ferromagnetism and SILC. Because < T < 1, gate im- 
balance (\pf—Pb\ > 0, so that sin 2 9 decreases and cos 2 9 
increases) makes the double-layer system more unstable 
towards pseudospin polarization. 

3. Subband densities versus bias 

In this section, we show some illustrative calculations 
of the effect of layer imbalance. We plot the subband den- 
sities n as and the value of L versus the gate imbalance 
parameter £ = (pf — PB)/n<T at fixed layer separation 
d/ao = 5 assuming zero interlayer tunneling (t = Q), for 
different values of r s . We find that if we fix layer density 



(or equivalently, r s ) and vary the gate imbalance param- 
eter C, then (for t = 0) there are six distinct patterns of 
transitions between noncrystalline MFA ground states. 
We list the six possibilities below, in order of increasing 
layer imbalance (beginning at £ = 0) from left to right, 
and in order of increasing average interparticle separation 
per layer r s from top to bottom: 

1. r s < r s (3,4): 

(p = 4)-»(p = 3)-»(p = 26), 
all with sin 9 = 

2. r s (3,4) <r s < r s (2,3): 
(p = 3) - (p = 2b), 
all with sin 9 = 

3. r s (2,3) <r s < V2r s (2,4): 

(p = 2o)-»(p = 3)-»( P = 26) J 
all with sin 9 = 



4. V2r s {2,4) <r s <r 
(p = 2o)-»(p = l) 1 
all with sin 9 = 



(1,2) and Cc > (1-Ti): 



5. %/2r s (2,4) <r s < r s (l,2) and Cc < (1 - Ii): 
(p = 2a) -> (p = 1), 

with sin6» > for Cc < C < (1 - ?i) 

6. r s > r s (l,2): 
(p = 1) only, 

with sin9 > for < C < (1 - Ii) 

Here (p = 2a) denotes the spin-polarized two- 
component state, (p = 26) denotes the pseudospin- 
polarized two-component state, and Cc is the value of 
the gate imbalance parameter C at which the (p = 2a) — > 
(p = 1) transition occurs. The quantity ^/2r s (2,4) ap- 
pearing in cases (4) and (5) above correspond to the crit- 
ical density n c for the MFA spin polarization transition 
for a single layer, expressed in terms of the average in- 
terparticle spacing per layer: l/\/ 7r ( ri c/2)a§. (Note that 
for the same total density ut, a double- layer system has 
an average layer r s = 1/ ' ^n{nT /2)a 2 that is \/2 larger 
than the single-layer r s = l/y^TrnxOg.) The quantity n c 
is discussed below Eq. (p2|). We shall illustrate the first 
four of these possibilities in the remainder of this sec- 
tion, and discuss the last two possibilities (which exhibit 
SILC) in Sec. VI. It is evident from the above list that 
SILC, which requires sin0 ^ 0, occurs only for p = 1. 

Figure [)] is an example of case (1), with r s = 1. This 
gives a four-component (p = 4) phase when the gates are 
balanced (C = 0), and maintains a p = 4 state for most 
of the range of C, followed by a p = 3 state for C near 
one. If we were to increase C beyond one (not shown), 
corresponding to pu < 0, then a pseudospin-polarized 
p = 2b state would be obtained. 
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FIG. 9. Normalized subband densities and interlayer ex- 
change parameter F, versus layer imbalance parameter 
C = {Pf — Pb)/tit, for d/ao — 5 and r s — 1. Fhe p — 4 
state is obtained except for £ near one, where the p — 3 state 
is obtained. Beyond f = 1 (not shown) the p — 2b state is 
obtained. 



Figure [Tl] is an example of case (3), with r s = 2.5. 
This gives a spin-polarized two-component (p — 2a) state 
when the gates are balanced (£ = 0), and maintains a 
p = 2a state for £ < 0.55, followed by a p — 3 state for 
0.55 < C < 0.95, and then a pscudospin-polarized p = 2b 
state for ( > 0.95. 
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Figure [lO] is an example of case (2), with r s = 
r s (2,4) m 2.011. This gives a three-component (p = 3) 
state when the gates are balanced (C = 0), and main- 
tains ap = 3 state for most of the range of followed by 
a pseudospin-polarized p = 2b state for £ near one and 
beyond. Note that for small £, T(p = 3) > 1. 



FIG. 11. Normalized subband densities and interlayer 
exchange parameter V, versus layer imbalance parameter 
C = {pf — pb)/tit, for d/ao = 5 and r s — 2.5. The p — 2a 
state is obtained for £ < 0.6, followed by a p = 3 state for 
0.6 < C < 0-8) an d then a pseudospin-polarized p = 2b state 
for C > 0.8. 




FIG. 10. Normalized subband densities and interlayer 
exchange parameter F, versus layer imbalance parameter 
C = (pf —PB)/n T , for d/ao = 5 and r s = r s (2,4) w 2.011. 
The p = 3 state is obtained except for £ near one, where the 
pseudospin-polarized p = 2 state is obtained. 



Figure [lj is an example of case (4), with r s = 9. 
This gives a spin-polarized two-component (p = 2a) state 
when the gates are balanced (£ = 0), and maintains 
a p = 2a state for £ < £ c ss 0.45, followed by a one- 
component (p = 1) state for £ > £ c . Because (c > 1 — T±, 
sin (9 = throughout, and thus no SILC is found. In 
the next section, we consider values of r s large enough 
that the p = 1 state is achieved for £c < 1 — Fi, thereby 
producing SILC. 
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FIG. 12. Normalized subband densities and interlayer 
exchange parameter V, versus layer imbalance parameter 
C — {pf — Pb)/tit, for d/ao = 5 and r s = 9. The p = 2a 
state is obtained for f < £ e » 0.45, followed by a p = 1 state 
(without SILC) for C > Cc- 



VI. PSEUDOSPIN-POLARIZED STATES 

In this section, we consider the case in which all the 
particles are in the lowest-energy subband [n a = fir), 
corresponding to full pseudospin polarization. There are 
two types of pseudospin-polarized MFA ground states: 
spin-unpolarized (p — 2b), and spin-polarized (p = 1). 
The spin-unpolarized case requires either interlayer tun- 
neling [t > 0) or gate imbalance \pf~ Pb\ > 0, or both. In 
the absence of tunneling, the pseudospin-polarized p = 2b 
state has sin 8 — (Ising-like pseudospin) , and occurs 
whenever the total density and the layer imbalance are 
sufficiently large. If the tunneling t is sufficiently large 
and the total density tit is not too small, then it is possi- 
ble in principle to obtain ap = 2b state with cos 9 = 0, for 
Pf = Pb- Tunneling also reduces the value of r s required 
to achieve the p = 1 state. Because the pseudospin- 
polarized p = 2b MFA state does not occur without bias 
or tunneling, we shall focus mainly on the spin-polarized 
one-component (p = 1) phase, which can in principle 
arise without bias or tunneling. The one-component 
phase is especially interesting because it can occur as a 
broken-symmetry ground state of a double-layer system 
in the absence of tunneling or layer imbalance, at very 
small particle densities and layer separations .113 

When n a = rip, the ground-state energy (p3|) becomes 



where T\ = T(n a = tit) and Q = (pf — Pb )/»t- The 
properties of Ti are described in Sec. A4 of the Ap- 
pendix. Equation (|9|) includes the two cases p = 1 
(n aT = n T ) and p = 2b (n aT = n al = n T /2). For t = 0, 
the p = 2b case has cos (9 = ±1, so we will focus on the 
spin- and pseudospin-polarized one-component (p = 1) 
ground state. 

A. Spontaneous interlayer coherence 

The pseudospin-polarized ground state offers the pos- 
sibility of SILC. Recall that SILC occurs when the off- 
diagonal (or interlayer) density matrix pi2 is nonzero in 
the absence of interlayer tunneling: 



P12 = y^(ciks c 2k«) = lj i ( n a- n b )sm6 ^ 0, 

ks 



(90) 



which requires both hnite pseudospin polarization (n a > 
rib) and sin# ^ 0. In the pseudospin-polarized ground 
state (n a = Ut), P12 is just the geometric mean of the 
layer densities: 



P12 = sin 6 = y/nin 2 , 



(91) 



where we have made use of Eq. (|T|) . So if the pseudospin- 
polarized ground state has some density of particles in 
each layer, it has interlayer phase coherence, p\2 7^ 0. 
Note that in the one-component phase, 



an f 



(P12 + P2i)/n T 



(92) 



measures the interlayer density matrix, normalized by 
the total density. 

For t = and p = 1, Eq. (f[S) gives 



m - n 2 



X = 



-1, X < —1 

= cos 6» = { X, -1 < X < 1 
I 1, X>1 

c C 



(93) 



1 - (32/3ttp)[1 - S(z)]/z 

c 

~^ (32/45tt)z - (l/24)z 2 ' 

where z = 2kpd, and we have made use of Eq. (53). 
The last line of Eq. (^) holds in the limit that z — * 0. 
The layer densities are equal {n\ = n 2 ) only at exacly 
Pf = Pb', when pf > Pb, layer 1 tends to be occu- 
pied, and when pp < Pb, layer 2 tends to be occupied. 
Thus, the hypothetical bistability of the one-component 
phaser-Hroposed by Ruden and Wu does not exist, due to 
SILciia Equation @ gives 
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sine = v /i-[c/(i-r 1 )] 2 e(i - r x - |ci)^ £ 



(94) 



so that sin (9 in nonzero only for |£| < 1 — Ti and for 
n a i = tit- It is interesting to note that when the ground 
state is pseuodspin-polarized, the dependence of the layer 
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layer imbalance £) 
of the electrons or 



densities on external parameters (e.g., 
does not involve the effective mass m* 
holes. 

We have found that layer imbalance {pp ^ ps) can 
induce SILC at higher total densities than in the balanced 
case. This is illustrated for r, = 11 in Fi g. |l3|, which is 
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Figure 
2a) state 



an example of case (5) introduced in Sec. 
|l3| shows a spin-polarized two-component (p 
when the gates are balanced (£ = 0), and maintains a p = 
2a state for £ < £ c « 0.28, followed by a one-component 
(p = 1) state for £ > £ c . Because £ c < 1 — Ii, sind > 
for £ c < £ < 1 — Ti, producing SILC in a finite region of 
layer imbalance away from £ = 0, at a smaller value of 
r s than is required to achieve SILC for balanced layers. 
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FIG. 13. Normalized subband densities, interlayer ex- 
change parameter F, and normalized interlayer density matrix 
sin (9, versus layer imbalance parameter £ = (pf — Pb) /tlt , 
for d/ao = 5 and r 3 = 11. The p — 2a state is obtained 
for £ < £ c ~ 0.45, followed by a p = 1 state with SILC for 
Cc < C < 1 ~ Fi. SILC is lost for C > 1 - IV 

When the toal density (and layer separation) are suffi- 
ciently small, the p — 1 phase with SILC is obtained even 
in the balanced case. This is illustrated for r s — 15 in 
which is an example of case (6) introduced in Sec. 
Figure |TJ shows a one-component (p = 1) state 
throughout the range of £, with sin# > for £ < 1 — Ti, 
producing SILC in that region. Note that as r s oo, 
then z — 2kpd — ► 0, so that Fi(z) — * 1, and thus 
1 — Ti — > 0. Therefore the maximum amount of im- 
balance ( which allows SILC decreases with the density, 
at very low densities. 
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FIG. 14. Normalized interlayer density matrix sm8 and in- 
terlayer exchange parameter V, versus layer imbalance pa- 
rameter ( — (pf — PB)jriT in the one-component phase 
(?iaT = n r), for d/ao = 5 and r s — 15. SILC is lost for 
C>l-Fi. 

Figure [l5| shows sin# versus C at d/ao ~ 4.356 for 
the three values of r s , where sm.6 1 obeys Eq. ([m|). At 
the lowest density shown (r s — 11), the system exhibits 
SILC when balanced (£ = 0) and under bias, until £ = 
1 — T\. As the density is raised, SILC is lost for the 
balanced system but appears suddenly around £ ~ 0.2 for 
r s = 10.5 when an abrupt intersubband charge transfer 
produces pseudospin polarization (p — 1). For r s — 9 
there is only a very small region of layer imbalance £ 
that exhibits SILC, and for r s slightly smaller than this 
value, SILC disappears completely. 
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FIG. 15. Interlayer exchange strength 

(pi2 + p2i)/riT = sin 9 as a function of the gate- imbalance 
parameter £ = (j>f — pb)/tit for fixed d/ao ps 4.356. 

Using Eqs. (|9|), (^ll), and (p3|), we may calculate the 
energy per unit area e\n\, n 2 ) defined in Eq. (^), for the 
one-component phase: 



£(711,712) = 



ILrp 



e 2 3/2 e 2 dni7t 2l 



47TC 



2c 



1) 



(95) 



where fif = 711+77,2. Recall that £ does not include the 
electrostatic contribution to the energy per unit area. 
We may use Eqs. (||) and (^5|) to calculate fj,i, the chem- 
ical potential measure relative to the energy minimum of 
layer 1, 
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where we have used the fact that in the one-component 
phase, Ti depends on 711 and 712 only through hit = 
n\ + ri2- The quantity «2 may be obtained by interchang- 
ing Tii and 712 in Eq. ( |96|) and ca n be used to compute the 
front-gate voltage Vf using Eq. (A3) in the Appendix. It 
is straightforward to check that the difference between 
fix and /12 satisfies the equilibrium condition in Eq. (|4|). 
Equations (§) and @ may be used to calculate the elec- 
tronic length Sjj defined in Sec. O: 
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The quantity S22 may be obtained by interchanging 711 
and 77 2 in Eq. @. 

Interlayer correlations produce a nonzero value of the 
electron length si 2 . From Eqs. (@), (0), and (p6h we have 



S12 



S » + 2 



Ti + (771 - 772) 



9nr 



Note that S21 = S12, and that 



si = sn - S12 



Tl + (771 - 77 2 ) 



g£i 

<9t7t 



(98) 
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The length S2 can be obtained from si by interchanging 
ni and 7i2- Note that 



si + S2 = -rid, 



(100) 



so that, from Eq. (|l l|) , the Eisenstein ratio for fixed total 
density (constant tit) is 



Re = 



-1 



(101) 



d + s 1-Ti (32/45tt)z' 

where z = 2fci?d, and the right-hand side holds in the 
limit z — > 0. It is interesting to note that i?g oc — 1/d 



as d — > 0, just as was found in Ref. |26| for the z/t = 1 
2LQH state. For fixed ps (nearly equivalent to keeping 
the back-gate voltage Vb constant), 



Re 



81 



[Ti + (77! - n 2 )dT 1 /dn T } 



(102) 



which in the balanced case (pf 
gives Eq. (|7q), which is exactly half of Eq. (101) 



2(1 -ro 

Pb so that ni = 772) 



Figure |16[ shows an example of the one-component 
phase under bias for fixed back-gate density pb (essen- 
tially fixed back-gate voltage Vb-) The normalized layer 
densities n\ and 712 are shown, together with the inter- 
layer density matrix pnjpB and T, and also the Eisen- 
stein ratio Re- For pf/pb < 0.5, layer 2 contains all the 
charge (772 = 77t) and layer 1 is empty, so that Eq. ( |i"o| ) 
gives Re = 1- For 0.5 < pf/pb < 2.25, both 771 and 
712 are partially occupied, P12 = ^Jn\n2 is nonzero, and 
Re has dropped abruptly and become negative, reflect- 
ing the prese nce o f SILC, with its value in this region 
given by Eq. (102). For pf/pb > 2.25, layer 2 is empty, 
77! = 777^, and Re = 0. Figure |l6| illustrates an interest- 
ing hypothetical situation in which bias and the exchange 
interaction have completely emptied out layer 2, despite 
the fact that ps is nonzero. It turns out that within the 
MFA, sufficently large pf would eventually repopulate 
layer 2. We analyze this issue further below, when we 
calculate the energy gap A a fc in a pseudospin-polarized 
state (n a — tit) to an otherwise empty subband [rib = 0). 
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FIG. 16. Layer densities n\ , 72,2 and normalized interlayer 
density matrix 2pi2/nr = sin#, in the one-component phase, 
versus front-gate density pf/pb, for fixed back-gate density 
Pb- 

Although Zheng and co-workers showed (within an un- 
restricted HFA) that an abrupt interlayer charge transfer 
does not occur when the gates are electrostatically bal- 
anced (pp — Pb), one may ask if an abrupt interlayer 
charge transfer can occur if the system is biased. Ac- 
cording to the MFA developed here, the answer is yes, 
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except at sufficiently small densities. (We are speak- 
ing here of zero tunneling; otherwise, neither layer is 
strictly empty, according to the MFA.) In the MFA, suf- 
ficiently strong bias or tunneling or sufficient lowering 
of the densities eventually produces an abrupt intersub- 
band charge transfer, i.e., at some point, n a suddenly 
jumps. Whether this translates into abrupt interlayer 
charge transfers depends on what happens with the phase 
angle 8. If sin^ 7^ (SILC), then the interlayer charge 
transfer is suppressed, or at least somewhat reduced. If 
the density is so high that strong bias produces the p = 2b 
pseudospin-polarized phase with sin 6 = 0, then the MFA 
does give an abrupt interlayer charge transfer, because 
then cos 6 = ±1, and there is no difference between sub- 
band densities and layer densities. So, for example, a 
system with a density that would correspond to p = 4 
when balanced will not exhibit SILC. It is likely that in- 
cluding correlation-energy effects eliminates the abrupt- 
ness of the transition, but these effects have not been 
included here. 

B. Intersubband gap 

The intersubband energy gap A Q b for the pseudospin- 
polarized (p — 1 or p — 2b) phase is defined as the energy 
required to move a particle from the occupied a subband 
with n a = nx to the (otherwise empty) b subband: 



for 



Aab = A (Jl_ 

Sn \ L X L, 



n a T ~> (l/p)(nr - 5n), Sn, 
n al -> (1 - l/p)(n T - Sn), n H = 0. 



(103) 



(104) 



An o utlin e of the MFA calculation of A Q b is given in 
Sec. A 6 of the Appendix. In units of the energy scale 
vq = e^/Aneao, the MFA intersubband gap is 
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(j>F - Pb)/tit, and 



2kpd is the layer 



where £ = 
imbalance. 

The intersubband gap A a b is useful for at least two 
purposes. First, it provides an estimate for the location 
of the pseudospin-polarization transition. The condition 
A a b > means the pseudospin-polarized ground state 
is stable against intersubband charge transfers, whereas 
A a b < implies the opposite. Thus, solving the equa- 
tion A ao = in the MFA yields an estimate of the loca- 
tion of the pseudospin-polarization transition. It turns 



out that this procedure gives a lower value of r s for the 
p = 1 transition than the GRPA estimate (obtained using 
the pseudospin Stoner enhancement factor /): at d = 0, 
the MFA A ab calculation (for t = ( = = cosfl = 0) 

gives ra°^(l,2) = 7r/\/2, compared to (1,2) = ir 
from the GRPA. Figure [ll] shows A O0 versus d/ao when 
t = ( = cos# = for r s = 1,2,3,6. It is evident 
that A Q0 > only for sufficiently large r s , and that it 
decreases with layer separation d/ao- Negative values 
of A ao indicate regions where the pseudospin-polarized 
state is not stable. 




FIG. 17. Normalized intersubband energy gap A a j, for 
transferring particles from the occupied lower-energy sub- 
band a to the (assumed empty) higher-energy subband b, for 
r s = 1,2,3,6, 10. 

The intersubband gap has a simple form when t = 
sin6 = 0: 



Aafc 



- + -(i-|CI) 

P ao 



2 4 1 



P 7T r s 



(106) 



Equation (|106| ) provides an estimate of when one of the 
layers cotains all the particles. It tells us that for suffi- 
cently high total density (small r s ), A Q b < and both 
layers must be occupied, provided that 1/2 + (d/a )(l — 
ICI) > (which includes £ = 1, corresponding to ps = 0.) 
This is because the kinetic energy (and the Coulomb en- 
ergy, for |f I < 1), which favors occupying both layers, 
dominates over the exchange energy, which favors occu- 
pying a single layer, at high er d ensities. For example, 
for p — 2 and £ = 1, Eq. (106) shows that both lay- 
ers will be occupied even when ps = 0, provided that 
r s < tt/2. On the other hand, for p = 1 and sin# = 
(i.e., ICI > 1 

vided that r s > (ir/y/2)\l - 

situation shown in Fig. n6l which shows that layer 2 has 
completely emptied for pf/pb > 2.25. At higher values 



Ti), only one layer will be occupied, pro- 
> (tt/V2)|1 + (d/o )(l - ICI)]- This is the 
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of pf/pb (e.g, r s < tt/2, not shown), layer 2 would no 
longer be empty. 

The second use of the MFA calculation of A a fc is as 
a rough estimate of the minimum energy (thermal, or 
from photons) required to excite particles from the oc- 
cupied to the unoccupied subband in the one-component 
state. Perhaps this could be detected with sensitive heat- 
capacity measurements or by measuring microwave ab- 
sorption. 

C. Coulomb drag 

One very interesting feature of the one-component 
state with SILC is that it is expected to exhibit inter- 
layer drag (finite dc transresistance) , even at zero tem- 
perature. Ordinarily, if the layers are not correlated in 
the ground state, current in one layer can drag along 
particles in the other layer (due to the Coulomb intex=. 
action between the layers) only at finite temperature.til 
But if pi2 ^ 0, due to either tunneling or, more inter- 
estingly, to SILC, then interlayer correlations present in 
the 2LES ground state will produce interlayer drag even 
at zero temperature, as has,been predicted for interlayer- 
correlated 2LQH states.LrEil Based on the Kubo formula 



with the ground state of Eq. (15), we expect that a cal- 
culation of the zero-temperature dc transconductivity ad 
will give 



e 2 pi2 e 2 (n a - n b ) sinfl 

(Td OC (X . 



(107) 



Calculations of the drag conductivity for a pseudospin- 
polarized ground g±3te are currently being carrie d ou t by 
other researchers. E3 According to Eqs. ( |9l| ) and (107), we 
expect that in the p = 1 phase (n a = tit), <7d ~ \/^i"2, 
approximately (i.e., within an MFA calculation of ad-) 

It would be interesting to clarify the relationship be- 
tween S12, P12, and ad- We conjecture that finite inter- 
layer drag at zero temperature requires pxi ^ at zero 
temperature (although we have not proved this) and that 
P\2 7^ (or at least P12P21 7^ 0) at zero implies finite in- 
terlayer drag at zero temperature. It is certainly true 
that P12 7^ and ad ^ occur together at zerg temper- 
ature in the interlayer correlated 2LQH effect .0 We also 
think it likely that a similar relation holds between S12 
and P12, and hence between S12 and ad, although we have 
not proved this either. 

VII. CONCLUSIONS 

We investigated the effects of intralayer and interlayer 
exchange in biased double-layer systems, in the absence 
of a magnetic field. This was accomplished using a mean- 
field approximation (MFA) which, in the limit of bal- 
anced layers (no bias), is equivalent to the unrestricted 
HFA of Zheng and co-workersED 



A. Findings 

We found that a balanced 2LES possesses four pos- 
sible noncrystalline MFA ground states. The spin- and 
pseudospin-unpolarized four-component (p = 4) state is 
obtained at the highest densities. In contrast to earlier 
work, we found that as the density is lowered, a three- 
component (p = 3) state with slightly unequal layer den- 
sities is obtained. Thus, we find that there is no direct 
four- to two-component transition. At finite layer sepa- 
ration (d > 0) and zero interlayer tunneling (i = 0), the 
(p = 4) — > (p = 3) MFA transition involves a small but 
abrupt interlayer charge transfer. Any such abrupt in- 
terlayer charge transfer will in principle result in a large 
(formally infinite) value of the Eisenstein ratio Re at the 
transition. The Eisenstein ratio is a sensitive measure of 
the interlayer capacitance discussed in Sec. H. 

Like Zheng and co-workers, we found that as the to- 
tal density was lowered, a spin-polarized two-component 
(p = 2a) state preceded a low-density one-component 
(p = 1) state possessing SILC, provided that the gates 
were balanced (pf = Pb)- This p — 1 state is differ- 
ent from that of Ruden and Wu, whose proposed one- 
component state occupied a single layer, rather than a 
single subband consisting of a linear combination of both 
layers. We obtained a MFA phase diagram for the non- 
crystalline phases of the 2LES, shown in Fig. |[ This 
phase diagram is similar to that of Ref. 17, except for the 
presence of the p = 3 phase between the p — 4 and p = 2 
phases. Only the p = 1 phase was found to possess SILC 
— i.e., a nonzero interlayer density matrix (p\2 ^ 0) even 
with zero interlayer tunneling (t = 0). We also defined 
the pseudospin Stoner interaction parameter J, and con- 
sidered the linear response of the MFA ground state to 
interlayer tunneling, equivalent to a GRPA calculation. 
We used / to obtain an alternate (GRPA) estimate of 
the location of the (p = 2) — > (p = 1) transition, shown 
as the dotted line at the top of Fig. a. Of course, in the 
limit of vanishing total density (r s — > 00), we expect that 
a Wigner crystal state is obtained (in the absence of dis- 
order). We did not consider the effects of disorder here, 
except to note that it limits the maximum r s for a state 
with mobile particles. 

Under bias (\pf — Pb\ > 0), we found that there are 
five possible noncrystalline ground states. In every case, 
the MFA gave subbands, which, when occupied, were ei- 
ther completely spin-unpolarized or fully spin-polarized. 
Including correlation-energy effects would likely produce 
ground states with intermediate spin polarisations, as 
is apparently the case in three dimensions .ell The ad- 
ditional state that can appear under sufficiently large 
bias and/or interlayer tunneling is a pseudospin-polarized 
two-component {n a ^ = n a | = tit/2) state, which we la- 
beled p = 2b. The p = 2b state requires bias and/or tun- 
neling, and has sin^ = (and thus no SILC) for t = 0. In 
Sec. ^ we studied the effect of bias by considering sytems 
at fixed total density nr, for a range of values of the 
layer- imbalance parameter C = (pf ~ Pb)/tit- We enu- 
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merated the six possible scenarios for bias-driven tran- 
sitions between the (five possible) noncrystalline MFA 
ground states for t = 0. We also showed that a very sim- 
ple model that assumes no interlayer exchange and no 
spin polarization is capable of fitting experimental SdH 
data quite well (see Fig. [(]) , and that a simple LDF model 
can do the same in the presence of interlayer tunneling 
(see Fig. |). 

We studied the one-component phase under applied 
bias, finding that bias lowers the r s required for SILC 
(see Fig. |l5). Within the MFA, SILC occurs only in the 
one-component phase: when t = 0, sin# is nonzero only 
when n a f = rip. Perhaps including correlation-energy ef- 
fects would allow SILC for states with partial pseudospin 
polarization. But p = 1 is only a necessary condition for 
SILC, not a sufficient one. We found that if the layer im- 
balance parameter was too large (Q > 1 — Fi), then SILC 
was lost. When SILC occurs, the MFA gave a value of the 
interlayer density matrix equal to the geometric mean of 
the layer densities: p\% = ^/nin 2 - For the case that SILC 
is present, we calculated the layer densities (rij), local 
values of the chemical potential (Hi), electronic lengths 
(sy ), and Eisenstein ratio (Re) (see Fig. fl6| ). 

Ruden and Wu originally predicted an abrupt inter- 
layer charge tranfer for t = at sufficiently low densities, 
and layer separations, in the balanced case (pf — Ps)-t3 
Like Zheng and co-workers, we found that an abrupt in- 
terlayer charge transfer does not occur in the balanced 
case, due to SILC. However, the MFA (and the unre- 
stricted HFA of Ref. [l?]) does produce abrupt intersub- 
band charge transfers, even for the balanced case. This 
a feature of the HFA that requires correlation-energy ef- 
fects to remedy. In the case of nonzero bias, intersubband 
transfers are equivalent to interlayer transfers if sin 6 = 0, 
which is the usual case, except possibly for p = 1. In- 
terlayer subband transfers at t = are reduced or sup- 
pressed in the MFA only for the (p = 2a) — > (p = 1) 
transition, and only if £ < (1 — Ti) in the p = 1 phase. 
So SILC, when present, does reduce or eliminate abrupt 
interlayer subband transfers, but the MFA does not al- 
ways eliminate them under bias (pp ^ Pb)- If the system 
is at sufficiently low density and layer separation that it 
stays in the p = 1 state, then there are no abrupt inter- 
layer charge transfers under bias in the MFA, despite the 
fact that the layers can empty out as pp is changed (see 
Fig. 13). 

We also calculated the intersubband gap A a j, for the 
pseudospin-polarized (p = 1 or p = 2b) phases within the 
MFA, defined as the energy to move a particle from the 
lower energy a subband to the higher energy (empty) b 
subband. This energy provides an estimate of the single- 
particle intersubband gap in the pseudospin-polarized 
p = 1 and p — 2b phases, and can be used to esti- 
mate the stability of those phases. If the p = 1 phase 
can be obtained experimentally, A a b might be measured 
using heat-capacity or microwave/optical techniques. A 
very interesting feature of the one-component phase with 
SILC is that it should have nonzero interlayer drag, even 



at zero temperature, with the size of the interlayer drag 
conductivity being proportional to the interlayer density 
matrix pi2- 

Pseudospin polarization can be detected by SdH mea- 
surements, which exhibit oscillations that are periodic in 
1/H (where H is the applied magnetic field). The peri- 
ods of the SdH oscillations are given by 
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where A as is the cross-sectional area of the Fermi surface 
perpendicular to the applied magnetic field for electrons 
in subband a with spin s. Knowing the total density 
np (e.g., from Hall measurements), SdH measurements 
of the subband densities n as could allow a determination 
of the degree of spin and pseudospin polarization. For 
example, in the case of equally balanced layers (n\ = 
712 = np/2) having a p-component ground state (p = 
1,2,4) in the absence of tunneling (t = 0), there is a 
single (p-fold degenerate) SdH oscillation period, 
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which allows p to be determined directly from SdH mea- 
surements. 

B. Can it exist? 

Can the one-component state be realized in the bal- 
anced case? The p = 1 state is a legitimate solution of 
in the HFA, but we have not examined its stability here. 
It is hypothetically possible that the p = 1 state might 
always be preempted by a Wigner crystal state. Conti 
and SenatoreEa carried out diffusion Monte Carlo (DMC) 
simulations in the d = limit, calculating the p = 4 
ground-state energy as a function of r s and using previ- 
ous single-layer DMC resultsEa to estimate the p = 2 and 
p = 1 energies. They also estimated the ground-state 
energy for a Wigner crystal state, and found that the 
p = 4 state is obtained for r s < 42, and that the Wigner 
crystal state is obtained for larger values of r s . In their 
calculation, neither the p = 2 nor the p = 1 states are 
ever favored energetically. 

Although the DMC results in Ref. [l8] show the need to 
examine the existence and stability of the p = 1 state be- 
yond the HFA, they do not rule out its existence. This is 
because at d = the fermions possess CP(3) symmetry 
(spin and pseudospin fully rotateable and interchange- 
able) , and estimating the p — 2 and p = 1 energies using 
single-layer results misses part of the correlation energy, 
which lowers the p = 1 and p = 2 ground-state energies. 
Ideally, a DMC simulation of CP (3) fermions would be 
most useful to determine theoretically if the p — 1 state 
can be obtained, but such calculations might be pro- 
hibitively difficult to carry out. As a start, allowing for 
the possibility of Wigner crystallization (broken trans- 
lational symmetry) within the MFA calculation would 
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be helpful. Alternatively, a time-dependent MFA calcu- 
lation of the collective mode would indicate where (for 
what density and layer separation) the collective mode 
of the 2LES goes soft (u> — > 0), signaling the onset of 
Wigner crystallization. We are currently developing such 
a calculation of the collective mode. Better yet would be 
a double-layer STLS calculation allowing for the possibil- 
ity of Wigner crystallization, or at least a determination 
of when the STLS collective mode goes soft in a double- 
layer system. 

More important is the question of whether SILC can 
be achieved experimentally in the absence of a strong 
magnetic field, which serves to quench the kinetic en- 
ergy of the particles. The MFA and HFA underestimate 
the value of r s required for the transitions, perhaps by 
a factor of 10. For example, the spin-polarization tran- 
sition foij— a single layer has been estimate to occur for 
r s ~ 20. Ea Although such high values of r s have been 
achieved in p-type GaAs samples, even higher values of 
r s will be required to achieve spontaneous pseudospin 
polarization. Disorder imposes further constraints, be- 
cause it limits the maximum r s for which the particles 
are still mobile. However, given the impressive progress 
in producing double-layer systems with ever-lower parti- 
cle densities and ever-higher mobilities, it does not seem 
prudent to rule out the possibility that such a state might 
someday be realized. _ . 

As pointed out bx-. Zheng and co-workersEZI and by 
Conti and Senatorejla most of the considerations pre- 
sented in Ref. (?] for SILC in the quantum Hall regime 
should be relevant to the p = 1 phase (with SILC) in 
the absence of a magnetic field. In both cases, there is a 
ground state with broken £7(1) symmetry (when t = 0), 
due to interlayer exchange at small layer separations and 
particle densities. It is therefore expected that the p = 1 
state will exhibit many of the novel features of the 2LQH 
state with SILC, including zero-temperature interlayer 
drag, vortex excitations [of the angle <j> in Eq. (|2^)] and an 
associated a finite-temperature Kosterlitz-Thouless tran- 
sition (for t = 0), and interesting many-body effec±s_in 
tilted magnetic fields for finite interlayer tunnclingjjia 

C. Speculations regarding the 2LQH regime 

It is interesting to speculate about the applicability 
of these ideas to the 2LQH effect at total filling factor 
unity.El For sufficiently small distances (d <C d c w 1.2£, 
where I = y/h/eB is the magnetic lengthQ), the 2LQH 
system exhibits a quantum Hall effect J3 Theoretically, 
the small- d 2LQH state has nonzero pi2 even for t = 
(SILG\_and the 2LQH system exhibits strong Coulomb 
drag.LrEl At sufficiently large layer separations, it is found 
experimentally that the quantum Hall effect disappears,!!!! 
and it is has been proposed that there, is a quantum phase 
transition to a state without SILC.Q The nature of the 
ground state for d > d c is a topic of active investigation. 
It has been analyzed as a system of two weakly coupled 



layers of v = 1/2 composite fermions (CF's). Theoreti- 
cal calculations of the drag at low temperatures predict 
that the drag resistivity should scale with temperature 
as T 4 / 3 , based on calculating the effects of gapge-Jluc- 
tuations on two CF layers in the metallic stateE3~E3. It 
has also been proposed that the weak coupling between 
the CF's in different layers produces BCS .pairing be- 
tween them at sufficiently low temperatures] 3 ^ and that 
this paired staff .leads to a finite drag resistivity at zero 
temperatureOcj 

We point out here that besides the apparent BCS insta- 
bility between CF's in different layers, double-layer CF 
systems might be unstable to pseudospin polarization. In 
the limit d — > oo, the double-layer vr — 1 system may 
be regarded as a p = 2a phase (i.e., spin-polarized but 
pseudospin-unpolarized) of CF's in zero effective mag- 
netic field [y = 1/2 per layer). Naively, the presumably 
large effective mass of the CF's would correspond to a 
much larger effective value of r s than for the zero-field 
case, perhaps producing a value of r s sufficently large 
to obtain a pseudospin-polarized p = 1 phase. Another 
way of saying this is that the large magnetic field ex- 
perienced in each v — 1/2 layer quenches the kinetic 
energy of the particles, and that this quenching might 
strongly enhance exchange instabilities - in this case to- 
wards pseudospin polarization (n a > rib), presumably 
with sin# = 1 (when pp — pb), since no spontaneous 
interlayer transfer has been found to occur for pp ~ pb- 
One appealing feature of a hypothetical p = 1 CF state 
is that it would have small or zero resistivity in the pseu- 
dospin channel (i.e., for oppositely directed currents), the 
same channel in which the 2LQH state exhibits superflu- 
idity at sufficiently small layer separations.!] Such a p = 1 
state of CF's would exhibit SILC (pu ^ even when 
t = 0) and therefore possess zero-temperature Coulomb 
drag. Recent experiments with double-layer systems cor- 
responding to filling factor v = 1/2 in each layer provide 
evidence for the possibility of zero-temperature drag £21 
We are currently investigating the possibility and con- 
sequences of pseudospin polarization in double-layer CF 
systems. 

We also note that a perpendicular magnetic field B 
will generally enhance spin and pseudospin polarization 
because it tends to quench the kinetic energy. This effect 
enhances the exchange and correlation effects that lead 
to polarization. We therefore expect that SILC can in 
principle be found at any filling factor vp — (h/e)riT/B, 
provided that the total density and the layer separa- 
tion d are sufficiently small. In particular, if it exists in 
the one-component phase for zero magnetic field, SILC 
will probably persist, and even grow stronger, when a 
perpendicular magnetic field is applied. Finally, we re- 
mark that it may prove instructive to view the vp = 1 
2LQH state as a Chern-Simons bosonic condensate of 
spin- and pseudospin-polarized (p — 1) electrons bound 
to unit flux quanta. 

It has recently been found that double-layer systems 
in strong magnetic fields near total filling factor unity 
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exhibit a huge resonant enhancement of the— interlayer 
tunneling conductivity when SILC is present .e9 It would 
be interesting to measure the tunneling conductivity 
for a tilted sample, since many-body effects in a state 
with SILC strongly suppress the interlayer tunneling am- 
plitude when the parallel component of the magnetic 
field exceeds a critical value. This suppression is much 
stronger than for a system without SILC. We expect a 
similar strong enhancement of the tunneling conductiv- 
ity at zero magnetic field, provided that the system pos- 
sesses SILC. Such tunneling measurements could prove 
very useful for measuring the strength of SILC in double- 
layer systems at zero (or higher) magnetic field. 
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APPENDIX: 

1. Gate voltages 

The layer densities (ni,Ti2) are determined theoreti- 
cally by minimizing the total ground-state energy per 
unit area [Eq. (p3|)] for fixed gate densities (pf,Pb)- Ex- 
perimentally however, it is the gate voltages which are 
tuned. If needed, the gate voltages for a given value of 
(pf,Pb) can be calculated using 



eVp = eDpEp + //1 + eV l 



(o) 

F > 



(Al) 



eVs = eD B E B + M2 + eVj 



(o) 

B ' 



and Eq. (g). Here the gate electric fields depend on the 
gate sheet densities through Gauss's law, E a = (e/e)p a , 

and eVa°^ are sample-dependent constant gate-voltage 
shift s. A lthough eV a is approximately equal to eD a E a , 
Eq. ( |Aip shows that the layer values ^ of the chemical 
potential also contribute to the gate voltages. 

The layer values of the chemical potential ^ can be 
computed numerically from the variation of the equilib- 
rium value of the total energy per unit area (regarded as a 
function of the front- and back-gate densities pf and ps) 
with respect to infinitesimal changes in gate densities: 



For the typical case in which the back-gate voltage Vb 
is kept constant and the back-gate distance Db is much 
larger than the interparticle and interlayer separation so 
that pb is nearly constant, it is convenient to use Eq. (^) 
and write 

eV F = eD F E F + edE 12 + A t 2 + eV F ° ] , (A3) 
where Eq. (|A5) gives 



d 
dp E 



Jo_ 

L X L V/ pr 



(A4) 



Equation ( [A3] ) has the advantage of being app licable even 
when layer 1 empties out. Equations (A3) and ( A4) 
can be used to calculate theoretically the front-gate volt- 
age. In the limit where the interlayer separation is larger 
than the intralayer particle separation and (which usu- 
ally amounts to the same thing) interlayer correlations 
can be neglected, then variations of \xi with pp are small 
in comparison with edEyi-, so that the effects of fi 2 can 



be absorbed into the voltage shift eV F , thus giving 



V F ^ E F Dp + E 12 d+V F 0) . 



(A5) 



2. Hartree-Fock approximation 

In the Hartree-Fock Approximation (HFA) the two- 
body interaction is factored so that the ground-state en- 
ergy per unit area is 



HF 



E X Ly L X Ly 



\~^2 £ k (4k s c lka + 4 ks C 2ks ) (A6) 



L x L t 



XMks^ks + 4k s C lk S 



jlkiSij2k 2 S2 
X ( C ]ik2Sl C J2k 2 S2)( C ]2klS 2 C Jl k l ; ''l) 

3l 32 

- X! X! V 3 a ( q = )P<* n 3 

j a 

+ 2^2 V "P(Q = °)P<xP0- 

a/3 

The effect last three (the Hartree) terms of Eq. ( |A6p may 
be calculated by noting that 



lim ViM = lim [l - (l - e'^)] (A7) 



S£ /L x L y = hiSpf + V2&PB- 



(A2) 



= lim da = (00) <L,-, 

\q^2eq) 2e 13 y 1 2e lJ ' 
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where (oo) d eno tes the formally divergent part in the last 
line of Eq. (|A7]). The last three (the Hartree) terms of 
Eq. (|Aq) become 



(oo)~(m + n 2 - pf -Pb) 2 + : tt 1 (pf - ni)(n 2 - pb) (A8) 



1 

'2 
e 2 D F 
2e 



PF^ni + n 2 - pb) + 



e 2 _d 
2c 
e 2 D B 
2e 



Pb{ui + n 2 - pf) 



Requiring the first term of Eq. (A8) to not diverge im- 
poses charge neutrality: ni+ n 2 — pf+Pb- From Gauss's 



law, the Hartree energy [Eq. (A8)] may therefore be writ- 
ten, up to an overall constant, as 



[E 2 2 d + E F D F + E 2 B D B ] , 



(A9) 



where E\ 2 is the electric field between layers 1 and 2, E F 
is the electric field between the front gate and layer 1, 
and Eb is the electric field between the back gate and 
layer 2. 

Equation ( |A9| ) is just the electric field energy per unit 
area for the sample; we drop the last two terms since they 
may be regarded as constants for fixed pf and ps- The 
ground-state energy per unit area may thus be written 
as 



£0 

L X Ly 



L X Ly 

t 



ks ci ks + c 2ks c 2ks ) 



(A10) 
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3. Exchange integrals 

The exchange integral I a p s (q) is defined as 



x6(^ s -|K-q/2|), 



L X Ly K 



where a and f3 can be either a or 6, and where k as , kp s 
denote the Fermi wave vectors for parti cles of spin s =| , | 
in subbands a or b. Note that Eq. (All) implies that I a /3 S 
is l/(27r) 2 times the shaded area shown in Fig. where 
for concreteness q is taken be in the x direction. 




FIG. 18. The quantity I a p s (q) is proportional to the 
area of the shaded region of overlap between two circu- 
lar Fermi surfaces of radii k aa and kp a , which are centered 
at K x — ±5/2. The circular Fermi surfaces intersect at 
K x = Ko/2. 

Let the quantity Ko/2 equal the value of K x at which 
the Fermi circles of radius k as and kp s intersect: 

K = (k 2 as - k 2 )/q = 4Tr(n as - n/3,)/q. (A12) 



Then 



Ial3s{q) = [nf3 s Q(k as - kp s - q) 

+ n as Q(kp s - k as - q)] 
+ Q(k as + k !js - q)Q(q - \k as - kp 3 \) 

1 r r -ifq + Ko 

x-{n as cos — 

7r V 2k ns 



(A13) 
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2k 
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q-K 
2k 



q-K 
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When P = a, then kp s — k as , K — 0, and Eq. (A13) 
becomes 



Iaas(q) = n as Q{2k as - q) 



(A14) 



COS 



2& n 



and the first exchange integral in Eq. ( p3[) may be carried 
out explicitly: 



1 8 
2L^E^(3)WS) = -^ 



,3/2 



Aire 



(A15) 



Equation ( A15 ) is just the exchange energy per unit area 
for a uniform spin-polarized two-dimensional electron gas 
of areal density n as . 
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4. Interlayer exchange parameter 

A key quantity in our discussion of the effects of inter- 
layer exchange in double-layer systems is the interlayer 
exchange parameter T, defined by 



r 



1 



L X Ly 



E 



Vll(q)-Vl2(q) 
e 2 d/2e 



(<?) + hbs(q) - 2i a bs(q) 



(A16) 



(n a - n b ) 2 



which is positive and a monotonically decreasing function 
of the interlayer separation d: 



< T(d > 0) < r(d-> 0), 



(A17) 



when the subband densities are regarded as fixed. We 
note that in the three-component phase (p = 3), unlike 
the (p = 1,2,4) phases, the equilibrium subband densi- 
ties change with d, so that according to Eq. (fig), 



lim T cq (p — 3) oc 



1 



(n a - n b ) 2 d a a 



oc > oo, 



(A18) 



whenpp = ps, in apparent disagreeme nt wi th Eq. (A17). 
We stress that the inequality in Eq. ( A17 ) is true only 
when the subband densities n as are regarded as fixed, 
which is not the case in Eq. (A18). 

The interlayer exchange parameter T is important be- 
cause it determines when SILC is possible (for £ < 1 — 
and p = 1) - i.e., when sin0 0. It also determines 
the value of the pseudospin Stoner enhancement factor I 
(which depends on To). T affects the state of the system 
(e.g., layer densities and Eisenstein ratio Re) whenever 
sin6> ^ 0. 



a. Inequality 

Using the inequality 

e 2 d/2e>V 11 (q)-V 12 (q), 



(A19) 



which is true for d > 0, it follows from Eq. ( A16| ) that 
for d > 0, 



r < 



L X Ly 



qs 

E 



Iaasjq) + hbsjq) - Zlabsjq) 

(n a - n b ) 2 
{n as -n bs ) 2 ^ r(d ^ Q); 



(A20) 



(n a - n b ) 

where we have used the fact that 

l ^[Vn(g)-VM) =1 
d^o e 2 d/2e 

and 



(A21) 



-j— ^2 [Iaas(q) + hbs(q) - 2I abs (q)] 



L X Ly 



(A22) 
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(^^j qKs 

x [Q(k as - |K - q/2|) - 6(fc bs - |K - q/2|)] 

E { 7^77 E t ( fc - - k ) - 6 ^ - w 1 

s K x y k J 

^(n as - n bs ) 2 . 



Thus the condition 

(»M - n b] ){n ai - n H ) > 



(A23) 



is sufficient t o gua rantee that T < 1. 

Equation (A22) is true at finite temperature when 
Iaf3s(q) is generalized appropriately. This is because 



li m [Vn(q)-V 12 (q)] = ^ 
d—>o ze 



(A24) 



is independent of the wave vector q. We may write 



1 



(L X Ly) 



E 

ki ,k2 ,S 



Vii(|k 2 -ki|)-yi2(|k 2 -ki| 

e 2 d/2e 



(A25) 



(< 



(&k 2 A 2 s>) (( a k lS a k lS ) - (&k ia &kl«) 



(n a - n b ) 2 

which generalizes V to finite temperatures. In the limit 
d — > 0, Eq. (A24) shows that T approaches 



1 



(n. Q 



E 



i 



(n Q - n b )' 



LjTy^H^ 



(A26) 



1, p = 1 (riot = «t) 

1, p = 2a (n aT = rt bT = n T /2) 

1, p = 3 (n aT = n ai = n fcT = n T /3) 

1/2, p = 2b (n aT = rt ai = n T /2) 

1/2, p = 4 (n aT = n Qi = n fcT = n H = n T /4) 

Empirically, we find that within the MFA, the spins 
in subbands a and b are either completely unpolarized 
(at higher densities) or fully polarized (at sufficiently low 
densities.) Therefore the only possible MFA configura- 
tions of spin and pseudospin that would not satisfy the in- 



equality in Eq. ( A23) and that might therefore have T > 1 



would be three-component (p = 3) states in which sub- 
band a (the majority subband) is spin-unpolarized and 
subband b (the minority subband) is spin-unpolarized: 



M = n a i,n a > > n bl = 0. 



(A27) 
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b. Pseudospin-unpolarized T 
It is useful to define and evaluate 



r 



lim r. 



(A28) 



We begin our calculation of To by noting that 
has (q) + hbs (?) - Habs (q) = (A29) 
J2 M k *s - l K + q/2|) - &(kbs - |K + q/2|)] 



ix-^J/ K 



x [9(fc as - |K - q/2|) - 6(fc 6s - |K - q/2| 
(^An/n r ) 2 



L X Ly 



■J2$(kF -\K + q/2\ 



K 



X5(fcjr - |K - q/2|) 
fc F An\ 2 0(1-1) 

27T 7lT / X\/l - X 2 ' 



where An = (n a — n ) — > 0, fcj? = ^jAwnr/p is the Fermi 
wave vector per layer for the state with p components 
(p = 2,4), and x = q/2kp. Using Eq. ( A29| ), we obtain 



2 2/tt 
J- o = - 



p e 2 d/2e J 
2 2 r^ 2 



^(W^ (A3Q) 



X" 



de- 



l 



suit 



where z = 2kpd, and the second line is obtained by the 
substitution of variables x — smQ. It is straightforward 
to obtain T p(z) for small z — > by expanding the last 
line of Eq. (A30) in powers of z, 



limr (z) = - 

2^0 p 



1 1 9 

1 - -z+—z 2 

7T 12 



(A31) 



up to second order in z. 

Obtaining To for z — > oo is more cumbersome. One 
way to proceed is to define a cutoff e that satisfies 



1 



< e < 1 



(A32) 



so that sin( 
Then 



r « — 
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P 7TZ 



for 9 < e and zsin6> 3> 1 for > e. 

e- ze ) 



(A33) 



r/2 , 

d6 



sin I 



The first integral in Eq. (A33) may be carried out using 
the identityBla 



f R (1 - p - *) 

/ dt [ - '- = ln(i?) + 7 + 



oo -t 

dt . 

t 



(A34) 



where 7 ~ 0.5772 is Euler's con stant , and R = ze — > 
00, so that the last term of Eq. (A34) can be dropped 
for large. R. The second integral in Eq. ( A33| ) is well 
known O 



d6- 



sm t 



ln[tan(6>/2)]. 



(A35) 



Now, tan(e/2) w e/2 for e < 1, so that the logarith- 
mica lly d ivergent (~ lne) parts of the two integrals in 
Eq. (|A33| ) cancel each other, leaving 



lim r = -— [ln(2z) +7] 

z^oo p 7TZ 



c. Pseudospin-polarized T 
We now compute 

Ti = lim T. 

n a —tllf 



(A36) 
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When the double-layer system is pseudospin polarized so 
that n a = t it, t hen I{, bs (q) = Iabs(q) = 0, and it follows 
from Eqs. (|A14|) and (A16|) that 



1 16 r , , 

Y\ = / dx (1 — e 

P nz J Q 

1 32 
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where z = 2kpd, and 
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1 - (3tt/32)z + (l/15)z 2 - (tt/256)z 3 , z -> 
37r/4z, z — > 00 



so that 







If 1 - (32/45tt)z+ (l/24)z 2 , z 
[Z) ^ P 1 (32/3^)/z-8/z 2 , z-^00 



(A40) 



Here Ji and Li are modified Bessel and modified Struve 



vtlff 



functions of the first kind, respectively:! 

I n (z)=r n J n (iz), L n {z)=i^ n+r >H n (iz), (A41) 

where J n is the ordinary Bessel function of order n, and 
H n is the Hankel function of order n. 



Obtaining Eq. (A3E) is somewhat involv ed. T he first 
part of the integral in the first line of Eq. ( A39 ) can be 
obtained by writing 



dxe zx arccos(a;) = / dxe zx arcsin(a;) 

Jo L 2 



(A42) 
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and using the identityB 

/ dxe~ zx arcsin(a;) 
Jo 



(A43) 
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[lo(z) - Lq(z) - e~ 



which corrects a misprint in Eq. (4.551.1) of Ref. 41. One 
then obtains 



/ dxe zx arccos(a;) 
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The second part of the integral in the first line of 
Eq. can be obtained by writing 



( dxe~ zx xVl - x 2 = ~ [ dxe- zx ^Jl-x 2 (A45) 
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and by using the identitiesil 
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to obtain 



dxe- zx xy/l - x2 = l + 7^ [L 2 (z) - h(z)} (A47) 



Combining. Eqs. (A44) and (A47) and using the 
identities^ 



I (z)-I 2 (z) = 
L (z) - L 2 (z) 



2J_(z) 

z 

2Lt(z) 2z 
z Sir 
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gives the second line of Eq. (A39). The third line of 
Eq. (A39) follows from power series (for small z) and. 
asymptotic (for large z) expansions of h{z) and Li(z).l__l 
We note that for the same value of p, Ti (z) < Tq(z) for 
all z. This may be s een by the subsitution of variables 
x = sin0 in Eq. flA3S| ): 



2 2 r /2 (l-e- zsin( 
Flit) = — / dd { - —— 



(A49) 



x {sin(26>) [(tt - 26») - sin(20)]} < T (z), 

where To is expre ssed as an integral over 9 in the last 
line of Eq. (A30). However, for differing values of p, 



Txip = 1) > T (p = 2) for < z < z c w 44.09. Both 
Ti(p = 1) and T (p = 2) are plotted in Fig. For 
z > z c (not shown), Ti(p = 1) < T (p = 2). 
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FIG. 19. = 1) (solid curve) and T (p = 2) (dashed 

curve) as a function of z = 2kpd. For < z < z c ~ 44.09, 

ri(p = i) >r (p = 2). 



5. Stoner enhancement factor 

We now outline our calculation of the Stoner in- 
teraction parameter for t — > 0. Consider a two- 
or four-component state with equal subband densi- 
ties (pseudospin unpolarized) , n a — rib = tit/2. 
For small t — > 0, imagine moving a small amount 
of charge An/2 from subband b to subband a, 
n a nx/2 + An/2, nt, riT/2 — An/2, so that (n a — 
rib) = An, and calculate the change in the total energy 
[Eq. (p^) ] as An — > 0. The effect of the change of den- 
sities on the last term of Eq. ( p3| ) may be calculating by 
using Eq. ( A30 ). The change in the energy per unit area 
due to An as An — > is then given by 



A£ (An) 



tAn 



1 



(An) 2 



pvo 
e 2 d f An 
~27 \ 2 



irp Aire 



(A50) 



r . 



Minimizing A£q/ L x L y with respect to An to solve for 
An and using the definition Eq. (|5^) of the Stoner en- 
hancement / gives 



2-Ktkp 
vo 



[\-p-k F dY^ (A51) 
[Vn{2k F x) - V 12 (2k F x 



2Vn(2k F ) 



dx- 



which is just the first line of Eq. (|5^). Equation ( |A51| ) 
is equal to the t — > limit of the Stoner interaction pa- 
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rameter calculated in the GRPA, given in Eq. (14) of 
Rcf. E23L 



where £ = (j>f — pb)/it-t and z = 2kpd. E quat ion ( A56| ) 
is expressed in dimensionless form in Eq. (105). 



6. Intersubband gap 

The subband transfer energy A ab for the pseudospin- 
polarized (n a = tit) p = 1 or p — 2b phase is defined 
as the MFA energy required to move a particle from the 
occupied a subband to the (otherwise empty) 6 subband: 



6 



£ 



OTL \^L x Ly 

where S£q denotes the change of £q under 

n a -i ^ (l/p)(n T - Sn), n 6T -> Sn, 
n al -> (1 - l/p)(n T - Sn), n H = 0, 



(A52) 



(A53) 



where Sn <C np, and n b i = reflects the fact that at 
low densities (n b = Sn — > 0), subband b will be spin- 
polarized. In order to calculate A ab we first compute 

Y^i 1 ^) + hbs{q) ~ 2I abs (q)} (A54) 
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x y K 



Sn L X L, 



Tn2 



x [Q(k as - |K - q/2|) - 6(fc 6s - |K - q/2|)] 
5] 6(fc 0fl - |K + q/2|) [Sk as 5(k as - |K - q/2[ 

Ks 
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on 27r z ' 

- n(Sk bs ) 2 e(k as - q)] 



-6(1 - x) arccos(x) + 29(1/2 - x) 
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s 1 



where fc as = -y/47m as so that 
<5fc as = -(2Tr/k F )Sn asl 



5k, 



'61 



V4vrfe, (A55) 



and fcf = \/ AnnT /p. In order to obtain the last line 
of Eq. (A.54), we expanded the bracketed terms to first 
order with respect to Sn, and defined x = q/2kp. 
The subband transfer gap is then given by 

A afc = - 

pp 
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